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Abstract 



We deal with the "nonrelativistic limit", i.e. the limit c — > oo, where 
c is the speed of light, of the nonlinear PDE system obtained by cou- 
pling the Dirac equation for a 4-spinor to the Maxwell equations for the 
self-consistent field created by the "moving charge" of the spinor. This 
^ ' limit, sometimes also called "Post-Newtonian" limit, yields a Schrodinger- 

^\ , Poisson system, where the spin and the magnetic field no longer appear. 

£~ — ■ However, our splitting of the 4-spinor into two 2-spinors preserves the 

' symmetry of "electrons" and "positrons" ; the latter obeying a Schrodinger 

m : equation with "negative mass" in the limit. We rigorously prove that in 

the nonrelativistic limit solutions of the Dirac-Maxwell system on R 1+3 
converge in the energy space C([0,T]; H 1 ) to solutions of a Schrodinger- 
Poisson system, under appropriate (convergence) conditions on the initial 
data. 

We also prove that the time interval of existence of local solutions of 

S Dirac-Maxwell is bounded from below by log(c). In fact, for this result 
we only require uniform H 1 bounds on the initial data, not convergence. 

Our key technique is "null form estimates", extending the work of 
Klainerman and Machedon and our previous work on the nonrelativistic 
' limit of the Klein-Gordon-Maxwell system. 

1 Introduction 



In this paper we study the behavior of solutions to the Dirac-Maxwell (abbr. 
DM) system in the limit c — > oo, where c is the speed of light. Coupled to the 
Coulomb gauge condition, this system has the form 

Wdp-M + rfA^p = o, d v F^ = j„/c, 9^ = 0. (i.i) 

AMS Subject Classification: 35Q40, 35L70. 
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Here the unknowns are the spinor field x/)(t,x) £ C 4 , regarded as a column 
vector, and the electromagnetic potential A^(t,x) € M, pk = 0, 1,2,3. Further, 
F^v = d^A v — dyA^ is the electromagnetic field tensor, and 

is the 4-current density. On the Minkowski spacetime R 1+3 we use relativistic 
coordinates x° = ct e K, x = (a; 1 , a; 2 , a: 3 ) G M 3 . <9 M stands for Thus, 
da = \d t , where d t = J^. We also write V = (<9i, <9 2 , <9 3 ), A = d\ + d\ + <9f 
and |V| S = (-A) s / 2 for sel Indices are raised and lowered using the metric 
(Vuv) = diag(— 1, 1, 1, 1). The Einstein summation convention is in effect. Thus, 
repeated greek indices /i, v, . . . are summed over 0, 1, 2, 3, and repeated roman 
indices j, k, . . . over 1, 2, 3. For example, A = djd^ . We denote by ( •, • ) c „ the 
standard inner product on C™. 

The physical constants are M — moc/h, g — e/hc, where mo is the spinor's 
rest mass, H is the Planck constant and e is the unit charge. By 7 M , fj, = 0, 1, 2, 3, 
we denote the 4x4 Dirac matrices, given in 2 x 2 block form by 

where the Pauli matrices a 3 are given by 

M^)- -"-(i-i)- 

The following related matrices occur frequently: 

where (fc, I, m) is any cyclic permutation of (1, 2, 3). Note the identities 

a j a k = -a k a 3 + 25^1 = S jk I + ie ]kl Si. (1.2) 

The first equation in Ijl.lfl is the Dirac equation. Multiplying it on the left 
by 7° and taking the imaginary part of its C inner product with ip yields the 
conservation law J M = 0. Thus, the "charge" is conserved: 

J (V>, VOc 4 dx = \\i>{t)\\ 2 L 2 — const. (1.3) 

The second equation in Ijl.lll is the Maxwell equation. We split A^ into its 
temporal part A , the electric potential, and its spatial part A = (A 1 , A 2 , A 3 ), 
the magnetic potential. Hence the electric field is given by E = VAo — 8qA 
and the magnetic field by B = V x A, and the second equation in 11.1(1 is seen 
to be equivalent to the Maxwell system in classical form, with charge density 
p = J°/c and current density ( J k )k=i,2,3- 
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The third equation in Ijl.lf) is the Coulomb gauge condition div A = 0. The 
reason for this choice of gauge condition will be explained later. It is equivalent 
to V A — A, where V is the projection onto divergence free vector fields in RjJ. 
The second and third equations in are then seen to be equivalent to 

AA - J°/c, {c- 2 d\ - A) A = c- Y V{ J fe ) fe= i, 2 ,3 

provided the initial data of A are divergence free. Thus, when properly rescaled 
(see pp, |201)) the system is conveniently expressed in terms of a small 

dimensionless parameter 

1 

£ ~ - 

C 

as follows: 

id t ^ = -ie^ojd^f + £r 2 7 V - A]a j ip £ - A e Q ^ E , (1.4a) 
AA e =p e , (1.4b) 
n E A e = eV3 e , (1.4c) 

where we have put in superscripts to emphasize the dependence on e. Here 

□ e = e 2 d 2 t - A 

and 

P s = (r,r) c *, J E = e - l {(«V e ,V' e > C 4} )fc=li2! 3- (i.5) 

We consider the Cauchy problem for (|1.4|l with "finite energy" initial data 

^ E \t=o = % e H\ (A £ ,d t A')\ t=a = (a§,a;) E VH l x Pi 2 . (1.6) 

We prove three types of results for this system as e — > 0. First, local well- 
posedness (abbr. l.w.p.) with a logarithmic lower bound on the existence time. 
Second, convergence in the nonrelativistic limit if the initial datum of con- 
verges. Third, we prove some more precise results on the asymptotic behavior 
of the Dirac spinor under various smallness assumptions on its "positron part" . 
These results are described in detail in the next three subsections. 



1.1 Local existence 

There are two issues here: (i) l.w.p. for e fixed, and (ii) the nature of the s- 
dependence of the local existence time as e — > 0. 

Concerning (i), the main difficulty is that one cannot directly estimate the 
bilinear term Aja J tp in the Dirac equation, due to the failure of the endpoint 
Strichartz estimate for the wave equation in 1 + 3 dimensions. The crucial 
fact proved here is that when the Dirac equation is squared, the bilinear terms 
resulting from this dangerous term can all be expressed in terms of null bilinear 
forms, provided the Coulomb gauge condition is used, and this enables us to 
prove l.w.p. of DM in the energy space (|1.6f> . a result entirely analogous to 
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that of Klainerman and Machedon f° r the Klein-Gordon-Maxwell (KGM) 
system. (The square of the Dirac eq. is similar to the Klein-Gordon eq., but 
contains some additional bilinear terms due to the presence of spin.) 

Bournaveas [S] proved l.w.p. of DM in the space (if>(t), A(i)) e H l / 2+s x 
H 1+d for S > 0, but this result does not take into account the null structure; 
in fact, by using the special structure of the equations and the so-called Wave- 
Sobolev spaces, the result can be improved to 2 (ip(t), A(t)) € H s x H 1 for 
1/4 < s < 1; this is proved in an upcoming paper by the third author. It 
is worth pointing out that these results are all independent w.r.t. s, since the 
regularity of A(t) is kept fixed. Thus, e.g., the l.w.p. in H 1 / 2 x H 1 does not 
imply l.w.p. in H 1 x H 1 (or vice versa). 

The question of global existence and uniqueness for DM remains largely 
open 3 (but see Georgiev j^j for a small data result), however, we prove — and 
this brings us to the second issue mentioned above — that as e — > the local 
existence time goes to infinity, subject to the initial assumptions 

1^611^=0(1), 11^11^+ e||af|| i2 = o( ±) as e — > 0, (1.7) 



s' 



where 



0<A<1 (1.8) 

will be kept fixed throughout the paper. (The upper bound 1/2 is explained by 
the factor e 1 / 2 appearing in the L 2 bilinear estimates discussed in Sect. [31) 

Theorem 1.1. (H 1 l.w.p. of DM.) The initial value problem \l-4\) , is 
locally well posed for fixed e, with an existence time T £ > depending only on e 
and the size of the norms of the data. Moreover, if l|1.7f) holds, then 

T s >c log- as e — ► 0, (1.9) 

£ 

where cq > is a universal constant, and we have 

U £ (t)\\ H i = 0(1), ||A £ (i)||^ + e ||d t A^)|| L2 = O (Jz) (1.10) 

uniformly in every finite time interval as e — ► 0. 

In order to control the evolution as e — > 0, it is crucial to have estimates 
which are sufficiently strong w.r.t. powers of e. To this end we employ analytical 
techniques used in our earlier paper [2], where the nonrelativistic limit of KGM 
was considered. The analysis of DM is more involved, however, due to the 
additional terms that come up in "squared Dirac" compared to the usual Klein- 
Gordon {KG) equation. In particular, we prove some new bilinear spacetime 
estimates which are needed to control these extra terms. 



2 Even this is not optimal (the scale invariant space is L 2 X H 1 ' 2 ), and in view of the recent 
work of Machedon and Sterbenz 1181 on KGM one may indeed hope to do better. 

3 This in contrast to the situation for KGM; see 1121 . The crucial point is that KGM has 
a positive Hamiltonian, unlike DM. 
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Once we have obtained closed estimates for the system — sufficiently strong 
w.r.t. powers of e — we use a bootstrap argument to prove existence in a short 
time interval depending only on the L 2 norm of iffi, provided e is sufficiently 
small, depending on the size of (|f .7|) . On account of the conservation of charge 
l|1.3|l for the Dirac equation we can then iterate this argument to obtain the 
long time result. 

We stress the fact that no convergence assumption is made on the data in 
the above theorem — all we need is the uniform bound l|1.7[) . However, if we do 
assume that i/'o converges in H 1 , then we can pass to the nonrelativistic limit, 
which we discuss next. 

1.2 Nonrelativistic limit 

The nonrelativistic limit of the linear Dirac equation with a given time-dependent 
electromagnetic potential was treated in pQ (earlier papers, see e.g. [H], dealt 
only with the static case, i.e. time- independent potential). There are also some 
results on the nonlinear Dirac and Klein-Gordon equations in the literature, 
see e.g. [221 ^ but f° r the coupled nonlinear Dirac-Maxwell and Klein-Gordon- 
Maxwell systems there are no results previous to our work (i.e. the present 
paper as well as [21 El) and the completely independent work of Masmoudi and 
Nakanishi [2T]. 

The most marked difference between our work and that of Masmoudi and 
Nakanishi is that our estimates are strong enough to give uniform (w.r.t. e) 
bounds for the solutions of DM assuming only the initial boundedness condition 
l|1.7[) — no convergence assumption is necessary. This, of course, is crucial as 
far as proving Theorems 11.11 and 11.71 is concerned. By contrast, in |22 the 
convergence assumption is essential because uniform estimates are obtained only 
on an arbitrarily small time interval, and in order to push the result to a larger 
time interval they must use bounds on solutions of the limiting Schrodinger- 
Poisson system. 

Let us now state our result. We split the Dirac spinor into its upper and 
lower components: 




(1.11) 



where x ancl V are 2-spinors, i.e. column vectors in C 2 . Before one can pass 
to the limit e — > 0, the rest energy must be subtracted, which for the upper 
"positive energy" component means multiplication by e lt / £ and for the lower 
"negative energy" component multiplication by e~ lt l e . 

Theorem 1.2. (Nonrelativistic limit of DM.) Consider the solution of 
{1.4V , if-? -6)) obtained in Theorem M.lV with data satisfying: 

(i) vo lim £ _,o 4*0 exists in H l , 

(ii) KH^+ellafH^ =0(£) ase^O. 
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Denote the upper and lower 2-spinors of vq by Vq and v respectively, and let 
(u, be the solution of the Schrddinger-Poisson system 41 



Au = n, n = \v+\ 2 + \v-\ 2 , ( id t ± — ) v± + uv ± = 0, (1.12) 



with initial data w±|t=o = v t ■ Then as e — ► 0, 

^ = e^( W +)+e^QQ+ (l) m H\ (1.13a) 

A £ =u + o(l) in H 1 , (1.13b) 

p e = n + o{\) in L p ,l<p<3, (1.13c) 

uniformly in every finite time interval. Moreover, the relativistic current density 
converges as follows: Let 

J° = Im ( Vv + ,v + ) C 2 - Im ( Vv_, u_ ) c2 

1 1 (1.14) 

+ -V x (av + ,v+) c2 - - V x (ctw_,w_) c2 

where (Vw±,w±) and (<rv±,v±) are the vectors with components (d : 'v±,v±) 
and ( a 3 v± , v± ), respectively, for j = 1, 2, 3. T/ien 

J £ — >J° m [^(Itxl^)]' weafc* (1.15) 

as e — ► 0. 

The first line in r.h.s. l|1.14|l is the conserved current associated to the limiting 
system (|1.12H , whereas the second line consists of the well-known divergence- free 
additional terms due to the interaction spin-magnetic field |17| . 

We can improve the convergence rate to 0(e) by strengthening the initial 
assumptions. Thus, we shall prove: 

Theorem 1.3. Strengthen the hypotheses of Theorem \l.<!\ by assuming 
11^11^=0(1), ||Va§|| H1 + e ||af|| H1 =of-^ 



and 



W «Q + 0(e) m H 1 



as e — > 0. Moreover, assume Vq G H 5 . Then 



V> £ = e~ lt/e2 ^+ j + 0(e) in H 1 as e — ^ (1.16) 

uniformly in every finite time interval. Furthermore, the convergence in Kl.lSb}) 
and \ 1.1 Sep is also 0(e). 



4 This system is globally well posed for L 2 data. 
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Remark 1.4. The hypotheses are not strong enough to guarantee strong con- 
vergence of the current density J £ locally uniformly in time. In fact, a simple 
counterexample is given by the initial datum 

* - G<) ■ 

Then J e initially has vector components 2 Re ( ct-^Vq, ) , which does not agree 
with the weak limit J° given by (|1.14|) . 

It is instructive to compare the last theorem to the formal derivation of the 
nonrelativistic limit usually reproduced in physics textbooks, the basic premise 
of which is a smallness assumption on the lower component ff of the spinor. 
The idea is to define 5 

«f = := e^V . (1.17) 

Then (|1.16Jl can be restated 

X e = v + + 0(e), rf = 0(e) in H 1 as e — > 0. (1.18) 
The Dirac equation i|1.4afl gives 

iD QX e = -iaWrf, iD oV e = -ia j D jX e - *rf, (1.19) 
where we write D = ed t — ieA e and Dj = dj — ieAj. Thus, 

rf = ~e\i^d jX e - e 2 \ {idtf + + Ay X e } , (1-20) 
and substituting this in the first equation in l|1.19|) gives, after some algebra, 
MtX e = \ (<V + eA £ ) 2 x £ ~ A%x e - \eB^x e ~ (1.21) 

where 

r e = l<T j Dj {d t if - iA&f) (1.22) 

and B e = V x A e . Then by formal considerations of magnitude, in particular 
assuming dtrf = 0(1), one obtains a Schrodinger equation in the limit e — > 0. It 
is possible to make this argument rigorous, but it has a fundamental weakness 
which limits its usefulness, namely that dtrj £ can be no better than 0(1 /e) unless 
one adds a further constraint on the initial data. In fact, it is clear from H1.20JI 
that d^f — 0(1) in L 2 initially if and only if the constraint 

?f = -e^ia j djX e + 0(e 2 ) (1.23) 

5 Here we break the symmetry of the signs in I1.13al . i.e. between "electrons" and 
"positrons", but this is not important since the lower component is in any case expected 
to vanish. 
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holds in L 2 at time t — 0, assuming the data Ijl.tifl are 0(1). 

However, the constraint l|1.23f) is not needed in Theorem 11.31 the reason 
being that instead of the simple splitting into upper and lower components as 
in (|l.llfl . we apply the eigenspace projections of the "free Dirac operator" 

Q £ = -iea k d k +-f°. 

As in pQ we use the spectral decomposition 

Q e = A £ n^ - A e ni 

where 

A e = \J \ — e 2 A, n| = ^(/±[A £ ]- 1 Q s ). (1.24) 

Since the positive and negative eigenvalues ±A e correspond to positive and 
negative energies of a free Dirac particle, the spectral decomposition is related 
to electrons and positrons ( 7 ). The formal limit e — > of TT§_ yields the 
operators 

n» ± = i(/ ± /), nl-(J S), n°_ = (° J). (1.25) 

The following basic lemma shows that is the leading order term in a series 
expansion of n|. in powers of e, and moreover that (|1.23|) is basically equivalent 
to H £ _ip £ — 0(e 2 ), a condition which resurfaces in the next subsection. 

Lemma 1.5. For all s € K, IIf|_ is bounded from H s — > _ff s uniformly in e. 
Moreover, 

n ±= n ± = F£^i (1-26) 

= nlTie^a k d k Te 2 n\ (1.27) 

where 1Z E - denotes an operator bounded from H s — > i/ 5- - 7 uniformly in e. 
Proof. This follows immediately from 

n £ ± - n° = T i[A e ]- Wft m- 1 ) 7° (1-28) 

and the fact that the Fourier symbol of 1 — [A 6 ]" 1 satisfies the inequalities 

0<1 <min{l,e|e|,£ 2 |e| 2 } (1-29) 

where £ is the Fourier variable corresponding to x. □ 

Before moving on, we prove that the initial data assumption 10) in Theorem 
1 1.21 implies something stronger, namely the convergence of IIj.i/jq. 
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Lemma 1.6. If 



in H 1 . 



lim = v 



exists in H , the 



limlT^o = ( (J J and lim Q u -' l l , o 



Proof. It suffices to prove (H§_ — H±) i/jq — > in ff 1 . But the proof of Lemma 
11.51 shows that the Fourier symbol of — LTj_ is bounded in absolute value 
by min{l,e|£|}. Thus (lT| - IT^) (-0g - v ) -> in if 1 , and by Plancherel's 
theorem and dominated convergence, (n|_ — II±) vq — » in If 1 . □ 



1.3 Semi-nonrelativistic limit 

As in pQ, by the "semi-nonrelativistic limit" we understand the approximation of 
the upper component of the Dirac equation by the Pauli equation for a 2-spinor, 
which reads 

idtXp = \ (<V + eA e ) 2 Xp - Af )X e P - \eB^ x p (1.30a) 
with initial condition 

XpI*=o = Xpq e J7 1 . (1.30b) 

Note that the naive "upper and lower components" approach in H1.19[) - I|1.22[l 
can give at best an 0(e) approximation to the Pauli equation, assuming the 
initial constraint l|1.23fl . which as remarked is essentially equivalent to H £ _iJ; £ = 
0(e 2 ). 

In contrast, by using the Dirac projections IIj_ instead of just ILj_, we can 
prove an 0(e 2 ) approximation, with the same initial constraint. In fact, we have 
the following result: 

Theorem 1.7. Consider the solution of \l-4\j , J^l.b}) obtained in Theorem M.lX 

Define \ e as in and let \p be the solution of the Pauli equation \1.3U\) . 

Assume the initial conditions 

(i) UlWm = 0(1), ||Va§|| H 4 + e ||af \\ H < = 0(1), 
(«) =0(e 2 ), 
as e — > 0. Then if 

llx £ "Xp|| ffl =0(e 2 ) (1.31) 

holds at time t = 0, it also holds uniformly in every finite time interval. For the 
current density we then have 

J £ = Jp + ^V x (a X p,Xp) C 2 +0(e) m L\, (1.32) 
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uniformly in every finite time interval, where 

3 P = lm((W-ieA') X £ P ,x P ) C 2 (1.33) 

is the current density of the Pauli equation. 

The remainder of this paper is organized as follows: In the following section 
we square the Dirac equation and reinterpret it in terms of the projections H^tp 6 
of the spinor, and we prove that the main bilinear terms can be expressed in 
terms of null forms. Then in Sect.|3]we discuss the linear and bilinear spacetime 
estimates of Strichartz type that are used in this paper. The proofs of those 
estimates that are not already in the literature can be found in Sect. ll2l In Sect. 

01 we define the function spaces that we use, and recall their main properties. 
The main estimates for the nonlinear terms are proved in Sect.El which is the 
heart of the paper. Then in Sects. I7HTT1 these estimates are applied to prove the 
main theorems. 

To close this section we introduce some notational conventions which will be 
in effect throughout: 

• For function spaces we use the following notation. If X is a Banach space 
of functions on K 3 , we denote by L\X the space with norm 

/DC \ l/p 

\\u(t,.)\\ p x dt) , 
-oo / 

with the usual modification if p — oo. The localization of this norm to a 
time slab St = [0,T] x R 3 is denoted \\u\\ L P X ^ ST y 

• In estimates, we use the notation < to mean < up to multiplication by a 
positive constant C independent of e. Moreover, in estimates over a time 
slab St, C is also understood to be independent of T. 

• For exponents, we use the convenient shorthand p + (resp. p~) for p + ( 
(resp. p—C) with £ > sufficiently small, independently of e. The notation 
oo~ stands for a sufficiently large, positive exponent. 

• We denote by f(x) <—> /(£) and u(t,x) i— > u(t, £) the Fourier transforms 
on M 3 and M 1+3 , respectively. As in [2] we split functions f(x) into their 
l° w (ICI ^5 an d high (|£| > 1/e) frequency parts, 

/ = /low + /high, (1.34) 

corresponding to a smooth partition of unity in Fourier space. 

2 Preliminaries 

As already mentioned, our approach to the Dirac equation is to square it and 
apply techniques similar to those used for KGM in . It is therefore convenient 
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to work with the "KG splitting" 

V4= (S) - l&^^T^idtr+AW)} (2.1) 

as used in 2 . In order to compare this to the Dirac projections 1)1. 24)1 . observe 
that if ip s solves the Dirac equation (|1.4a|l , then 

r± = ^ =p ^[XT'iA^r)- (2.2) 
But using the estimate 

||[AT7|U ^ £ ~ T ll/llfl-r for 0<r<l, (2.3) 
followed by Holder's inequality and Sobolev embedding, we see that 

£ 2 ||[A1-H^W)|| H . ^^"IIAII^II^II^ for 0<<7<1, (2.4) 

so r.h.s.JOJ is 0(e 1_A ) in i? 1 at time i if the bound lfl~TU|) in Theorem H~T1 
holds. As far as proving Theorem 11.21 is concerned, it is therefore immaterial 
whether we use ip± or Tl±ip E . 

For later use we note the following consequences of (|2.2|l and (|2.4|) . First, 

\H\\ H „<\m\ H *+e 2 -n\^ s \\m\m\m < a < 1, (2.5) 

using the uniform boundedness of n§_. Second, 

\\X-\\ Ll + \W+\\ Ll < e \W\\m + W A£ Wm \\1> e \\m > ( 2 - 6 ) 
where we used l)1.2fi)l and the orthogonality of IT+ and U°_ . 

Let us now restate the system (|1.4|l in terms of the splitting (|2.1|) of the spinor. 
First we subtract the rest energy, defining 



_ (X± \ ._ p ±it/e 



V4- (2-7) 



Thus 

^ = ^ + + 1 pt= e- lt ' e2 4> £ + + e +lt / e2 (t) £ _. (2.8) 

Lemma 2.1. In terms of the splitting 1)2. 8[l . defined via l|2.1)l and 1)2.7)1 . the 

Dirac equation {l.^a^ is equivalent to a system of two equations 

L%^ + = -A%<t>% + \^R S , L*_ 01 = -J§4L - ^e-^R*, (2.9) 

provided the constraint 1)2.2)1 is satisfied at time t — 0, or equivalently that the 
Dirac equation is satisfied at t = 0. Here 

L%=id t T^^ (2.10) 
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and R E is given by 

X s R E = e {2iA e • V + idivA £ + iE £ a j - B e S J }ii £ 

+ e 2 (A £ ) 2 ^ £ - [Ac, X £ W + - $L). (2.H) 

Further, [•, •] denotes the commutator and 

E = (E\E 2 ,E 3 ) := VA - ed t A, B = (B 1 , B 2 , B 3 ) := V x A. (2.12) 

Proof. Squaring the Dirac equation (|1.4a|l yields (cf. [7| Sect. 70]) 

{e 2 (id t + A £ ) 2 + (V - ieA £ ) 2 - e~ 2 - ieE^oP + eB £ S j } i) £ = 0. (2.13) 

Applying idt + A £ to both sides of 1(2.1(1 and making use of (|2.13() and 

e 2 {id t +Ai)r = x%r+-r-i 

which follows from (|2.1|) . one easily obtains (|2.9|l . Reversing these steps, one 
finds that l|2.9|) implies the squared Dirac equation 1(2.13(1 . But the latter implies 
the Dirac equation, since we assume that (|2.2|l holds initially, which amounts 
to saying that the Dirac equation is satisfied initially. □ 

Let us make a brief, heuristic comparison of 12.9fl with the expected limit 
11.12JI . As it turns out, R £ vanishes in the limit, so 1(2.9(1 tends to the Schrodinger 
equation in 1(1.12(1 . In fact, the Fourier symbol of (X £ — l)/e 2 is 

MO- ~('f/ 2 Tlf^r < 2 -»» 

,2 |£| e for |£| > 1/e, 



so L± tends to the Schrodinger operator idt ± A/2 as e — > 0. Moreover, the 
charge and current densities 1)1.5(1 are given in terms of the fields 1(2.7(1 by 

P E = ( X%, X% ) + ( Xt, Xt ) + < rf + , V e +) + ( if- , rjL ) (2.15) 
+ 2 Re {e- 2 ^ 2 < X E +) X I > + e"^ < r£, < >} , 

J £ = H Re {( ^ , ^ ) + ( < ) (2.16) 

+ e 2 "/^ 2 < ) + e- 2lt/£2 < ^X £ +,rf- >} . =1 2 g . 

We expect [cf. (I2.6(l ] that X-,?7+ — * 0. Thus, in r.h.s. ((2.15(l only the first and 
fourth terms are of importance, and A^4q = p £ tends to the Poisson equation in 

For later use we note the estimate 

0<|5|/e-M5)<e- 2 . (2.17) 
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This reduces to r — a(r) < 1, where 



Q(r) : = TTTTTt 1 ' <2 ' 18) 



But r — a(r) = r+ l — s=l where s = vl^ 



r 



2 



We now turn to the problem of obtaining closed estimates for the modified DM 
system (|2~3j) . l(Llb) . (JQcJ. A serious obstacle to estimating the bilinear terms 
in (|2.11|) is the failure of the endpoint Strichartz estimate for the wave equation 
in dimension 1 + 3. The salient feature of the Coulomb gauge, however, is that 
these problematic terms can be expressed in terms of the null bilinear forms 

Qo(u, v) = doudov — Vu • Vu, Q a p(u, v) = d a udpv — dpud a v, (2.19) 

where do denotes ed t and < a < (3 < 3. These bilinear forms enjoy better 
regularity properties than generic products of derivatives. 

We emphasize that in the following result ip docs not have to solve the Dirac 
equation. 

Lemma 2.2. (Null structure.) Given a potential {A^(t, x)} satisfying the 
Coulomb condition div A = 0, let E and B be defined as in (gJil), and consider 
the bilinear operator 

rp — > {2iA • V + iEjO? - B,jS J ) tj) 
appearing in \2.11)) . We have the following identities: 

2A-VV> = -Q J -fc(|VrV fc ,V) (2-20) 

and 



{;;/•;, i) r \;:\<v* n.s^c 

= QAM' 1 £dta jk , U) - QjkflVr 1 dia jk ,a l U) 
+ Q {A s ,cPU) + Q 0] {A k ,a^a k U) - ~Q jk (A mi e^ kl S ia m U) 



(2.21) 



Oik = RjA k - R k Ai, Rj = I VI 1 dj 



wher 

and U = f/(V>) is the J^-spinor defined by 

n e U = -i(s8 t +d>d j )il>, U|*=o = 0, isd t U\t= = Vb- (2.22) 
Here ipg denotes ip\t=o- 
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Proof. The identity l|2.2U|l goes back to the work of Klainerman and Machedon 
[T2"] on KGM, so we concentrate on the new identity (|2.21() . Define 

d± = ed t ± a^dj 

and observe that 

{d-Aj)a j = -(Ej - tf,.l,,:n' - ,7*,.S'\ (2.23) 

where we used the second identity in (|1.2|) and the assumption divA = 0. By 
the first identity in i|1.2fl . 

d+d- = D e . 

Thus Ij2.22(l implies that w = ip — id-U satisfies d+w = with w\t=o = 0, 
whence 

ip = id-U. 

Apply H2.23[l to this and use 

d + {a?U) = a j d-U + 2d j U (by JOJ) 
to rewrite l.h.s. l|2.21[) as 

{d-Aj)d+{ a iU) - 2{d-A J )d 1 U. 

To the last term we apply the identity l|2.20|l : this we can do since div9 M A = 
for fi = 0, 1, 2, 3 by assumption. To the first term we apply the following general 
formula, obtained using the second identity in (11.21) . 

(d-<j>)(d+U) = Qo(^U)+Q 03 (^a^U) - % -Q jk {(j>,e jkl SiU), 

where </) is a function and U a 4-spinor. This last formula is due to Klainerman 
and Machedon |H]. □ 

3 Bilinear spacetime estimates 

The main technical tools used in this paper are spacetime estimates of Strichartz 
type for solutions of the free initial value problems 

a £ u = 0, it| t=0 = /, d t u\ t=0 = 0, 
L±v = 0, v\ t=Q =g, 

on K 1+3 . Let us first describe the new L? product estimates that are proved in 
this paper, and then we recall the estimates proved in 0. 

Let fj, and A be dyadic numbers of the form 2 J , j 6 Z. Denote by A M the 
Littlewood-Paley operator given by 

(A M /H0 = /?(£/m)/(0, 

where (3 is a bump function supported in |£| ~ 1 such that X^ez Z^/^"') = 1 
for £ 7^ 0. We write = A p / and similarly for g, u, v. Thus / = etc 
We shall prove the following: 
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Theorem 3.1. The solutions u, v of 1)) satisfy the following dyadic spacetime 
estimates: 

(i) \\^{u x vx)\\ Llx < e 1/2 M Wfxhi \\9xL* if y. < A < l/e. 

(ii) \\A^u x v x )\y ^e^WWfxW^hxW^ if A » l/e. 

t.x 

(Hi) \\uu,Vx\\ L 2^ < e 1/2 min(/i, A) ||/ M || i2 \\gx\\ L 2 for all fi, A. 

See [SJ Thm. 12.1] for the analogous estimates in the case where u and v 
both solve the wave equation. 

By decomposing the product uv into dyadic pieces, then applying Theorem 
13. II and finally exploiting the orthogonality properties in Fourier space to sum 
up, one obtains the following corollary. (The complete argument can be found 
in [i Sect. 12].) 

Corollary 3.2. The solutions u,v of ^.1}) satisfy 

|||Vr (H|| L2 <C.x,. a e 1/2 ||/ll*. 1 \\9\\h,* 

t,x 

provided that 

si,s 2 <l, o < -, si + s 2 + er = l. 

Estimates of this type for the case where u and v both solve the free 
wave equation were first investigated by Klainerman and Machedon. The case 
(si, S2, cr) = (0, 1,0) is excluded, a fact related to the false endpoint case of the 
Strichartz estimates for the wave equation in 1 + 3 dimensions. However, by 
assuming a little extra regularity one can easily sum the dyadic pieces and one 
obtains the following nonsharp bilinear estimate. 

Corollary 3.3. The solutions u,v of \3.1\) satisfy 

\\uv\\ L ^<Css l ' 2 \\f\\ L2 \\g\\ H1+s 

for all 5 > 0. 

Proof. It suffices to prove the sharp estimate 

II^aII^ <e 1/2 \\f\\ L2 X\\gx\\ L2 . (3.2) 
Write u — ^2 and consider the cases fj, < A and /i>A. In the first case, 

where we used Theorem 13. If iii'l to get the last inequality. In the second case we 
have, by orthogonality in Fourier space, 

E, Ua] V\ < ^ ||UuUA||r2 j 

M>A ^ I ri Z-^/ n ~ i> \ >* "Jj 



L 2 ^m>A 



and by Theorem 13. lf iii) we dominate this by e ||/|| 2 2 A 2 ||<7a||l2- □ 
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Here we could also take / in H 1+s and g <E L 2 , but we shall not need this. 
However, for null bilinear forms one can get the sharp result (i.e. 6 = 0). Thus, 
we recall the following, proved in [2J Proposition 4]: 

\\Qij(u,v)\\ LU <£ 1/2 \\f\\H*\\9\\m (3-3) 

where Qij is given by (|2.19[) . We remark that this is the analogoue of an estimate 
for two solutions of the free wave equation proved by Klainerman and Machedon. 

Since we will prove part (|nj of Theorem 13. II by a reduction to linear Strichartz 
estimates, let us recall these (for 1 + 3 dimensions). We say that a pair (q, r) of 
Lebesgue exponents is wave admissible if (q,r) ^ (2,oo) and l/q+ 1/r < 1/2, 
and sharp wave admissible if the last inequality is an equality. 
For the free wave u in (|3.1|) one has the well-known estimate 

hW^^C^e 1 ^ (3.4) 

for wave admissible (q, r) and s — 3/2 — 3/r — 1/q. As proved in |16| . this can 
be improved if the Fourier support of / is small. Thus, if / is supported in a 
cube with side length ~ fi and at distance ~ A from the origin, where fi <C A, 
then 

\\u\\ LtL r<C q , r s^(^y /2 ' 1/r \\f\\ As , (3.5) 

for q, r, s as above. 

For v satisfying (|3.1|l we have, as proved in :2, Proposition 1], 

IMIifi; < Cq,r (l|flW||ffi/<> + ^ llffhigh || ^2/, ) (3.6) 

for sharp wave admissible (q, r). (Then one can use Sobolev embedding to obtain 
estimates for all wave admissible pairs.) In order to prove Theorem lOjujl we 
need the analogue of (|3.5|l in this context. Thus, we shall prove: 

Proposition 3.4. Let v be as in ^ff.ij) . and suppose g is supported in a cube 
with side length ~ fi and at distance ~ A from the origin, where fi <C A. Then 

\\v\\ LiL r < C q ,. (^) V2 1A (11^11^/,+ e^llflhighllfl*/.) (3-7) 

for sharp wave admissible (q, r) . 

Finally, recalling the basic heuristic that L £ ± behaves like a Schrodinger 
operator at low frequencies, it is not surprising that we have the following 
Schrodinger type estimates, proved in |2J. We say that a pair (q, r) is Schrodinger 
admissible if q, r > 2 and 2/q + 3/r = 3/2. 

Proposition 3.5. Let (q,r) and (q,r) be any two Schrodinger admissible pairs. 
Then for the solution of L±v = F with data u|t=o — f we have 

ll u low|| i 9 i r( ST) + \\v\ow\\ L oo L 2($ T ) < ||/low|| L 2 + \\Flovr\\ L f L r' (s^ ' 

where i + A- = 1 and i + X = 1. 

an r r 



16 



4 Function spaces 

We shall use the following spaces of functions on M 1+3 with weighted norms 
defined in Fourier space: 



• H s / with norm || (£) s (|t| - e~ x \£\} 8 u(t,0\ 

• H s e - e with norm |||£| s (|t| - e" 1 \£,\) 9 u(t, £)\ 

• Tip with norm + e\\dtu\\ H s-i,s. 

• Hl' e with norm ||u||^-i,e + e \\dtu\\ H o,e . 

• Kt±h c (i) with norm \\(0 s (rTh e (0) e u(T-,0\\ L 2 and h £ as in ISm 

Here (•) stands for 1 + |-|. These spaces are by now standard, and we will recall 
their main properties without proofs. For more details and further references 
to the literature, the reader may consult e.g. |23], [T5]- 
It will be convenient to introduce the notation 

U e (t) = e^-lVe 2 = e ithe(M)^ 

S(t) = e ItA / 2 , (4.1) 

W s {t) = e u W/ e 

for the propagators associated to, respectively, the operators L± defined in 
Lemma \2. II the Schrodinger operator and the wave operator. 



(i) Superposition principle. A fundamental property of the so-called "Wave 
Sobolev space" Hp 6 is that any function in this space can be written as a 
superposition (iP-valued integral over the real line) of solutions of the free 
wave equation with initial data in H s . (See Proposition 3.4] for the precise 
statement.) This, in effect, replaces Duhamel's principle in the framework of 
the Wave Sobolev spaces, and it has the following simple but extremely useful 
consequence (see [T5]h 

Transfer Principle. Suppose T is a multilinear operator (f±(x), . . . , fk(x)) 
T(/i, . . . , fk){x) acting in x-space. If T satisfies an estimate 

\\T(W*(±t)h, . . . , W%±t)f k ) \\ L , Li < Ce 1 ^ UMI^ ■ • ■ \\f k \\ H . k , 

for all combinations of signs, then 

||T (ui, ...,Uk)\\ L * Ll < C e e 1/q Huill^i,* • • • \\u k \\ H s k ,o 

holds for all Uj £ Hp' 6 ' , provided 9 > 1/2. Moreover, the same statement holds 
with H s and Hp ' replaced by their homogeneous counterparts. 
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The spaces X"'_ ±h are related to the equation L±v = in the same way 
that the Wave Sobolev spaces are related to the free wave equation. Thus, we 
have a superposition principle and hence a transfer principle for these spaces as 
well. To be precise, in the above Transfer Principle, one can replace any one of 
the W e (±t) by U e (±t) and correspondingly Hi 3 ' 6 by X^if^ Applying this 
to estimates from the previous section, we have, for 9 > 1/2, 

||'«high|| i g L r < e 1 / 2 + 1 / T " HMhighll^i.f for sharp wave adm. (q,r), (4.2) 

\\u\ovi\\ l i L t < ||uiow|| x o,e for Schrodinger adm. (q,r), (4.3) 

* ^ T = ±h e (£) 

||Wlow|| i 2 i cc < ||ltlow|| X M ■ (4.4) 

Here l|4.3l) follows from Proposition 13.51 with F = 0. By Sobolev embedding we 
reduce (|4.4H to the case (q,r) = (2,6) of (|4.3|) . Finally, l|4.2[) holds by virtue of 
l|3.6[) and the trivial estimate 

11/highlU <e'||/wBhlU.+- for ff>0 - ( 4 - 5 ) 

(ii) Embeddings. The most basic embeddings are 

H s /, X s r i ±he{0 ^ C b (M; tf s ), C* fc (R; ff s ), (4.6) 

which hold uniformly in e for any 8 > 1/2. Also uniform in e are 

L\Ll^ H°/-\ X° T £±l i() for |^<P<2, ~ < 6 < 1. (4.7) 

In fact, the dual statement H®' 1 ' 6 , -X^7_ A ^ «-> follows by interpolation 

between the trivial case p = 2, = 1 and (|4.6(l . We shall also need 

\\e ±u ^ 2 u\\ x o, e ^ < s- 2 ^ \\ u \\ . (4.8) 

"^x = ±f le (£) 

This is obvious if u(r, £) is supported in | |t| — |£ | /e| < e~ 2 ; then we can in fact 
replace the left hand side by ||w|| i2 - On the other hand, if u(t, £) is supported 
in |I T I - l£l l £ \ > e_2 5 tnen P^|> follows from l(2~T7|) . 

(hi) Time cut-off. In view of l|4.6|l . we can localize to any finite time slab 



S T = [0,T] 



x 



The restriction space H^' 6 (St) is complete when equipped with the norm 

II u IIhi' 8 (s t ) := inf {II^IIhi' 6 : v = u on ^} ■ ( 4 - 9 ) 
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Norms on the other restriction spaces are similarly defined. When 9 < 1/2 the 
embeddings (|4.6() fail, but since Hp 6 etc. are spaces of tempered distributions, 
it still makes sense to restrict them to the interior of St, and we will use the 
same notation Hp 9 (St) etc. for these spaces. Taking the inf over all extensions 
produces a seminorm in this case. 

The idea behind the following "cut-off lemmas" originates in the work of 
Bourgain 0] on the Schrodinger and KdV equations, and was developed further 
by Kenig-Ponce-Vega ^D] m their work on KdV and by Klainerman-Machedon 
[T5] and the last author [23 for the wave equation. In fact, the argument given 
in applies to X s ' 6 spaces in general, and in particular proves the following. 

Lemma 4.1. Suppose L±v = F on the interior of St with v\ t= Q = /. Let 
9 > 1/2. Then for < T < 1, 

ll«ll^ (e) (5,)^^(ll/ll fl . + Mx S - ft) (5 T) ) 

where C'e is independent of T and e. 

By rescaling x — > ex we reduce the next result to the case e = 1, which in 
turn follows from estimates proved in jl3J . 

Lemma 4.2. Suppose D £ u = F on the interior of St with (u, dtu)\t=o — (f,g). 
Let 9 > 1/2. Then for < T < 1 , 

\\u\\^e {ST) <C e {\\f\\m+^h\\^ + - e \\F\\ H y-\s T ^ ( 4 - 10 ) 

where C'e is independent of T and e. Also, for MeH large enough, 

II«II w ;.» (St ) < Cee- M (\\f\\ Rs + \\g\\ H .- x + \\F\\ nl -i,e-i {ST) ) . (4.11) 

The last inequality is not sharp w.r.t. e, but it will only be used in a situation 
where powers of e are not important. In order to estimate the Dirac current 
density we shall need the following "integration by parts" -version of Lemma 14.21 

Lemma 4.3. Suppose D E u = e %t l E F on the interior of St with vanishing data. 
Let 1/2 < 9 < 1. Then 



i|u|I 7^' 9 (St) ~ £ \\ d t F hHs T ) + \\ F ^t\\ L 2 H i +e\\(dt) F c: 



l^ext|| L 2 



for all < T < 1 and all extensions F cy ± of F to all o/IR 1+3 . Here (dt) is the 
multiplier with Fourier symbol (1 + 

Proof. Let us denote F cxt simply by F. Write 

e l */ £2 F = {e 2 /i)d t \e lt/e2 f] - (e 2 /i)e lt / £2 d t F 
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and u = i*i + u 2 accordingly. By i|4.1U|) . Huall^i.e^) ^I^-hI^^v Now 

define G = e < £ F. Split G — G\ + G 2 by a partition of unity in Fourier space 
such that 

Gi(t, is supported in |t| < |£| /e, 
Gi{T,g) is supported in |t| > |£| /e, 

and write u\ = iti,i + accordingly. That is, DeUij = (s 2 /i)dtGj on the 
interior of St with vanishing data. By (|4.10|) . ||iii,i|| 7 ^i,e, s \ < £ ||<9 t Gi|| i2 , but 
using Plancherel's theorem and the assumptions on the Fourier support, 

WWiL' z ^IH v i g i|Il 2 ^ ^IH v i g IIl 2 = JllM^IU 

whence Hui.iH^i.e^-, < II-^IIl 2 ^!- Finally, to estimate u± t 2 we first observe 
that it has an extension to all of R 1+3 defined in Fourier space by 

^(r,0 = _ £2 J +lel [(s 2 ^)dtG 2 ]-(r,0- 
Thus 121,2(7-, 01 ~ ]7| |G 2 (t-, 0|, and since 

hi,2\\ul- e {s T ) - Iki^ll^.e^i+s) < ||(£|r| + |f|)(|r| - |Cl/e) 8 ui,2(7-,0|| i 2 4 
we conclude that ||mi,2||^i>«(s t ) is dominated by 



e\\ (r) 9 G 2 (r, < e\\(r + l/e 2 ) 6 F(r, < s \\ (d t ) e F\\ J2 + » F[ ' 



This ends the proof of the lemma. □ 

5 Main estimates 

Here we prove the main a priori estimates for the nonlinear terms in the modified 
DM system, in terms of the following spacetime norms. 



Definition 5.1. For T > we define 

Iw'' 8 (St)' 



• X£ = e A ||A £ i 



Y T — Z)± P*± x 1,8 (St) 



• Z T — J2± (^±)low|| r a r ,8r 



lL?L8 nLrL 2 (ST) , 

for > 1/2 sufficiently close to 1/2, independently of e, but depending on the 
fixed parameter A. In fact, the relevant condition is 

A + l-2(9>0, (5.1) 

which we assume from now on. 
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We also need the following initial data norms. 
Definition 5.2. For initial data ll.Bfl we set 

. ^ = e A (K|| Al+£ ||af|| i2 ), 

• i? = WII HI , 

In order to simplify the notation we drop the superscript e on the fields 
<fi, ip, etc. in the remainder of this section. We assume < T < 1 in the 
estimates that follow, and we write 

P| = P(Xf + Yj) 

where P(x) = x + x N for a sufficiently large N £ N, independent of e. 

5.1 Estimates for A 

Split tp = V'low + V'high and write 

A = A' + A% (5.2) 
where A' corresponds to "low-low" interactions: 

AA' = (V'low^low) ■ 

Then 



AA' 


LfLL 3/2)+ (s T ) ~ ^ 2 


for 


1 <P< 2, 


(5.3) 


A< 


'l?l< 3 / 2 ' + (s t ) ~ £l ^r) 


for 


1 <p< 2, 


(5.4) 


AA' 


L 2 L^(S T ) ~ (-^t) 2 


for 


3 

!<*■<-, 


(5.5) 


AA'J 




for 


3 

l<r<-, 


(5.6) 


AA 


L 2 t L-(s T ) ~ (-^r) 2 


for 


2 < r < 6, 


(5.7) 


AA 




for 


1 < r < 3. 


(5.8) 



Here ()5.8JI follows from Holder's inequality and Sobolev embedding, while (|5.7|l 
reduces to 

\m LtLi r {ST) < Y£ for 2<r<6. 

By Sobolev embedding and the Transfer Principle, the latter reduces to the 
L\h\ Strichartz estimate (|3.6|) . Let us now prove 1)5. 3|) and 15. 4|) ; the proofs of 
<|5.5[l and (|5.6|l are similar. Write 

ll(V^)ll L (3/ 2)+ < m ht < mil w5 • ( 5 - 9 ) 
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For -0 = ip\ ow the L\ norm of this is clearly dominated by r.h.s. l|5.3[) . On the 
other hand, if at least one V'high is present, then we dominate by r.h.s. l|5.4[l using 
the H 1 <^-> U% Sobolev embedding and the estimate 14.5fl . 
We will also need the embeddings 



l/lli-. £ l|A/|| i( 3 /2) - + ||A/|| i( , /a)+ , 



IV/II 



< 



|A/|| 



|A/|| 



(5.10) 



5.2 Estimates for the remainder term 

For the remainder term R e given by (|2.11[1 we shall prove (cf. Lemma \A.l\i 



ix 



x = ±h e (5) (ST) 



< (1/2-A)- p e 



(5.11) 
(5.12) 



Using (|2.11[> . (|2.3|l . (|5.2() and 14.7(1 we dominate l.h.s. l|5.11|l by a sum of terms 



Ni = 

N 2 = 
N 3 = 

Ni = 

N 5 = 
N 6 = 



3="*/= - floaty - S,5^}|| x o^ k)(St) , 



!|A ■ V^|| i2(Sr) , 
\\e(d 3 A )a 3 ^\\ L 2 H1{STr 



-\\[A' Q ,X^ -l]i; ± \ 

£ " 1 

i||K,A E -l]^ ± | 

H £ ( a )vil 2(St) . 



L 2 (St)' 

L 2 t -Li(s T y 



All these terms appear also in the KGM case (see j2J), with the notable exception 
of N\. The latter is however the most interesting (and difficult) term, so we 
consider it first. Write {i(Ej — djA^a 3 — BjS^} tp = ^2 where 



{i(Ej - djA )a j - BjS j } 



and the sum is over all dyadic numbers fi of the form 2 J , j € Z. Here A^ t is the 
Littlewood-Paley operator defined in Sect. El We split into the cases 

(i) At < 1/e, 

(ii) ii > 1/e. 

Case 0). By ()4.7|l . we can reduce to proving 

< (l/2-A)- X ey 



M <l/e 



T i T , 



(5.13) 



£ 2 (S T ) 

but this follows from Corollary 13 . 31 via the Transfer Principle. 
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Case Using (|4.7|) we write 

||e ±It/e ^L«^ (0 (5 T ) ~ 11^11^(5.) \\ e±lt ' e ^WxHiI^St) ( 5 ' 14 ) 

where < a <C 1 will be chosen later. Proceeding as in case (0), but using the 
sharp estimate (|3.2[) . we obtain 

\\I,\\ms T )^ £l/2 ~ Ax TYS. (5.15) 
We claim there exist ( > and M G N, both independent of e and /i, such that 

h^lJ^-x <e" A V C ^- (5-16) 

Granting this for the moment, we see that by choosing a sufficiently small in 
H5.14f> . depending on M, we get 

E||„±rf/e 2 j II < ^.(1/2-A)" pe 

fi>l/e 

as desired. Let us prove the claim. On account of Lemma \'2. 21 
h = QjfcdVr 1 ed t a jk , AJJ) - Q^dVp 1 V*, cM^C/) 

+ Q (^-, a J 'A M C/) + Q 0j (A k ,a>a k AJJ) - % -Q jk (A m , e^S^A^U) 

where a jfe , C/ are as in Lemma 12.21 But since ip solves the Dirac equation, 

U e U = -e~ 1 'y ip + EA j oe i ip + eA ip (5.17) 
Now we appeal to the following null form estimate. 
Theorem 5.3. Let 1/2 < 6 < 1. Then 

\\Q{u,v)\\ H a.o-i < Cg IMI^i.e 

/io/rfs on R 1+3 for all null forms Q in \2. Moreover, if Q = Qij, then the 
norm ||u||^i,e in the r.h.s. can be replaced by \\u\\^i,e. 

By a standard procedure we reduce this to well-known bilinear estimates for 
the homogeneous wave equation; the proof can be found in Sect. H^l 

Applying this estimate, and recalling l|4.8|l . we reduce (|5.16|) to proving 

IIA^U^c,^ < C c e- M ^P^ (5.18) 

for £ > such that 1 + 8 < 2 — 2£. Clearly, it suffices to show 
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but using H4.11fl and l|5.17|) we reduce this to 

\\ A i aj *P\\ L *m-<(s T ) ^ Cce^X^YS, (5.19) 
PoV-|| L?ffl( 5 T) < (Yt) 3 - (5.20) 

The former follows from Corollary 13.21 and the Transfer Principle, while the 
latter reduces to (|5.8|) using Leibniz' rule, Holder's inequality and (|5.1U|) . This 
concludes the estimate for Ni. 

It remains to estimate the terms iVjj, . . . , Nq. Use Lemma f2. 21 and <|3.3[1 via the 
Transfer Principle to see that 

^2 < e 1/2 - K X £ T Y^. 
Next, by Leibniz' rule, Holder's inequality, l|5.10|l and l|5.7|l . 

To the term N4 we apply we apply the commutator estimate 

|| [A -1 (/S), A £ - l]h\\ L2 < e 2 C p \\f\\ m+P \\g\\ m+P \\h\\ Ht (for all p > 0) 
proved in j2J Lemma 9] . Thus 

iV 4 <e 1_ (y^) 3 . 

In N5 we simply expand the commutator and apply the estimate 

||(A £ -l)/|| i2 <e 11/11^, (5.21) 
which follows from I2.14|l . Thus 

N 5 < (||^'|| Lriso(Sr) + \\VAi'\\ LrLUST) ) \\4>±\\ LrsHST) , 

so in view of l|5.10|l , l|5.4|l and l|5.6[) , N5 satisfies the same bound as N4 . Finally, 
by Holder's inequality and the H 1 <—* L% Sobolev embedding, 

This concludes the proof of l|5.11|l . 

Now consider the estimate (|5.12|l . Using (|2.3|) with r = + we get 

WWl'(St) (Nl+N 2 +N 4 + N 5 +N 6 ) +N 3 , 

where N%,N3, N4 and N§ are as before, whereas 

Ni = \\i(Ej - d 3 A Q )aH - BjS^\\ L3tB0 - (ST) , 

^ = i||K,A s -i]^ ± || L?H0 - (ST) . 
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Write N\ < Ni^+Ni^ corresponding to ip — ip lam +V'hi g h- For the low frequency 
case we apply the nonsharp bilinear Strichartz estimate in Corollarv l3 . 31 via the 
Transfer Principle, to get 

Nx,i < s {1/2 - Kr X £ T Y£. 
By Sobolev embedding and Holder's inequality, 

Ni,2 < (l|VA|| £f , £ j +e\\d t A\\ L?L ^ ||Vhigh|| i? + L? - , 

where the pair (2 + ,oo~) is chosen to be sharp wave admissible. Applying the 
Strichartz estimate (|4.2|l we then obtain the same estimate for N\p as for N\,\. 

Next, 

^5 < (KII l? z~- (St) + ||V^'|| l?L 3 (St) ) U±\\ LrHHST) < e(YSf, 
where we used <|5.6[) to get the last inequality. This ends the proof of (|5.12|) . 

5.3 Estimates for the current density 

Split i\> — V'low + '/'high and write 

J = J' + J" 

where J' corresponds to "low-low" interactions: 

eJ' = {(a fe ^iow,^ow)} fc=1 ^ 3 . 

By |OJ| and 

||whighWlow|| L 2 < 1 1 "high 1 1 L oo L2 l^lowllijioo < e ||u high || x i,e ||viow|| x M , 

||whighVhigh|| L 2 ^ ||whigh|| L 4 ||^high|| L 4 < || Whighllx 1 ' 6 H^high || ' 

T = ±h € {£) T = ±h £ (£) 

whence 

I|J"IIl 2 (5 t) ^(^) 2 - (5-22) 
In order to estimate J' we expand it as in (|2.1(i|) . Thus, we write 

J' = (J')i + (J')a 

where 

e(J')i - 2Ref e - 2it / £2 < ^( X+ ) low , fa-Jw)} 

whereas (J')2 consists of products containing at least one of the fields (x-)low 
or (ry+)iow, which we expect to be small. The latter we estimate, just to take 
one of these terms, 

\\(v :l (X+)l0W, (tyf)loW )|| L 2( St ) < II(X+)1ow|| L 2 L oc (St) ||(??+)low|| L oo i 2( 5T ) • 
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Thus, using ll4"4"|) and 

\\m\ LHST) <(YS) 2 + e l - A P^. (5.23) 

To that part of A which corresponds to (J')i we are going to apply Lemma FOl 
Hence we want to estimate 

.. a ||F cxt || ^2 

Mt ■■= e ||a t F|| i2(Sr) + ||Pext|| iiffI + e ||<3 t ) 9 F cxt || i?:c + (5.24) 
where F ext is an extension of 

F = {(ct j (x+)iow,(j7-)iow)} j=12 3 (5.25) 
from St to all of R 1+3 . To choose this extension, let 

A - (£) e A-f ±Mf) 

be arbitrary extensions of <fi± and define F cxt by l|5.25|l with \+ & n d V- replaced 
by their respective extensions. From now on we denote F ext simply by F. We 
claim that 

M T <(Y T ) 2 +e {1/2r P T +e 1 ~ 29 U' + \\ xl . e UL\\ x i,» (5.26) 

If this holds, then taking the inf over all extensions yields 

M £ T < e l - 20 P T . (5.27) 

Let us prove (|5.26(l . First, 

£ \\ d tF\\ L 2 {ST) < - ||(A £ - 1)x+\\ ltl1(St ) \\(V-)iow\\ L *L~(s T ) 

+ £ \\ L +X+\\ L 2 (St) ll(»7-)lowlli~(S T ) + (•■■) ( 5 ' 28 ) 

< (F T ) 2 + e (1/2r {Y T f + £ (3/2r P|, 
where (...) stands for symmetric terms. Here we used (I5.21f) . <|4.4ll and 

\\L% X+ \\ L2(ST) <{Y £ f+eP T , 

ll(^-)low|lioo(5 T ) ^ £ ' 1 ^ 2 ' ) Yj- 

The former was obtained from l|2.9|l . (|5.12|l and (|5.8I) . while the latter follows 
from H4.6f) and Sobolev embedding. Second, we write 



i f Hl?//i < llxVILrn 1 IK 7 ?-) 1 ™ 1 1 + IKx'+Wll^oo \W- 



, II ' II I ' II ( 5 - 29 ) 



1,9 ! 



26 



where we used (14 .4f) and (|4.6|) . This estimate can of course also be used for the 
term ||F|| L 2 . Third, 



e||(a t ) e F ext |L 2 ^s^H^ILm Ik'-llvM , (5-30) 
where we used the following: 

Lemma 5.4. ||(d t ) e (ui ow vi ow )|| 2 < e~ e \\u low \\ i, e \\v low \\ i,e 

Proof. To simplify the notation, let us write u,v instead of wiow^iow here. 
W.l.o.g. we assume u(r, £),v(t, £) > 0. Then using Plancherel's theorem, 

\\(d t ) 9 {uv)\\ L2 < \\uT%\\ L2 +e' e \\u\V\ e v\\ L2 + \\vT + u\\ L2 + e- e \\v |V|%|| L2 . 

Here T± is the multiplier with symbol (r =F h e (£,)) e , and we used 12.14JI . Write 

WuT^vW 2 < |M| £f > £ s ||r*u|| L?L9 < \\u\\ x i,e \\v\\ x i,e 

using 14. 6|) and Sobolev embedding. Next, 

|h|V| e w|| 2 < \\u\\ LrL 3 |||V|%|| L2i6 < \\u\\ x i,o \\v\\ x i.» 

where (|4.3|l with (q, r) — (2, 6) was used. □ 
Finally, combining and (|53U|) . we get (|^|) . 

6 Iteration scheme and local existence 

For fixed e we shall prove the following local existence theorem: 

Theorem 6.1. For fixed e, the Dirac-Maxwell-Coulomb system (|1.4fl is locally 
well posed for initial data in the space (|1.6() . The existence time T > only 
depends on e and the size of the norms of the data, and the solution is in the 
space 

reHl' e (S T ), A^HY{S T ), A^eC^nH 1 ), (6.1) 

for all 1/2 < 8 < 1. Moreover, the solution is unique in this regularity class, 
and we have 

where <f>± is defined by 1|2.1[) and <|2.7|) . 

We shall prove this by Picard iteration. In order to simplify the notation we 
drop the superscript e on the fields tj), etc. and introduce instead a superscript 
(m) to denote the m-th iterate of a field. For (|1 -4ft we use the iteration scheme 

{ied t + icJdj - (l/e)7 } ^ (m+1) - -eA { ™ ] a 3 ^ {m) - eA { ™ } ^ m \ (6.2a) 

A4" l) =P M , (6.2b) 

□ eA (m+l) =e pj(m) j ( 6 .2c) 
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with initial data as in (|1.6f) . where p( m ) and j( m ) are given by 1)1. 5|) with tp 
replaced by its iterate ip( m \ Note that Aq is not really iterated; l|6.2bjl simply 
defines in terms of Observe also that for all to, 

divA^ =0, 

since w = A' m ' — PA'" 1 ' satisfies Dw — with vanishing initial data. 

By convention we start the iteration at m — — 1 and set all iterates identically 
equal to zero there. Then the iterates are just solutions of the free 

Dirac and wave equations with data (|1.6p . Define (cf. I|2.1II and H2.7[) ~) 

^ m+ l) = 1 |^ (ro+ l) ± £ 2 [A£] -1 ^ira+l) + ^l^Wj | ^ (6 3a) 

(6.3b) 



(m) •- c V± 



Proceeding as in the proof of Lemma 12. II one finds 



where 



S (m) = {2iA^ • V + i £ J (m W - B^ m) S j \ ^ {m \ (6.5) 
C (m) = {^4'™- Va fe + Af^A^-^aP - A^A^a^ V (m-1) , 

and E^ m) ,B^ n) are given by (|2~T2|) with ^ replaced by Alr^ 

We now turn to the proof of Theorem 16. II By standard arguments, this reduces 
to proving closed estimates for the iterates in the space (|6.1|l . Set 

rO) _ |L/,(m) II , || a (m) II 

and denote by Bq the norm of the data (|1.6|l . Then it suffices to prove 

B { ™ +1) < CP(B ) + CT S P (4 m) + 4 m_1) ) - (6-6) 

for some constants C, S > and a polynomial P with P(0) = 0. Here C and P 
may depend on e, but since the latter is fixed we do not indicate this explicitly. 
In what follows, C, 8 and P may change from line to line. (Observe also that 
since all the nonlinear terms in DM are in fact multilinear, the same arguments 
then give estimates for a difference of two iterates.) By Lemma f4. 21 

ll A(m+1) L- (ST) < CB ° + CT V m) \\l H s T y ( 6 - 7 ) 
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where the T comes from Holder's inequality in time. Now apply the Strichartz 
estimate ((3.611 via the Transfer Principle to see that ||V' (m) || i 4 (ST) < B^. In 
order to estimate -0( m+1 ) we use the splitting 1(2. 8|) and the embedding 

\\u\\ Hi ,e<e- 29 \\u\\ x .,o (6.8) 



^,(m+l) II .„ . < C /J|</' \\y'<- ,■ 

\\H E (St) Z-^M r =t nX r=±h e (o ( - ST > 



r = ±h c (0 

which holds in view of 1(2.17(1 . Thus, we write 

',("+!) II , „ <-r>^IU( m+li 
± 

Using Lemma ETT1 and (14.71). we bound ||0± llyi.e /„ by 

CP(B ) + CT s \\A^U i r ) \\ LrHHST) +C\\e^ (6.9) 
The second term is trivial to bound, since 

ll A 4 m) |L ri , (ST) <lk (m) ||^ (Sr) for l<r<3 (6.10) 

and 

||v4 m) (*)IL < ll^ m) (*)Lx + e ||A(™- i )(t)||^ 1 ||v (m - 1) (*)Lx- (6-ii) 

The latter is just the analogue of 12.5(1 for the iterates. 

For the third term in l|6.9|) we can apply the estimates proved in Sect. 15.21 
In fact, we claim that the proof of (|5.11() gives 

U e ±it/e* R (m)U <CT S p(B { T m} +B { ™- 1) ) . 

To see this, consider one by one the terms N%, .. . ,Nq in Sect. 15.21 For N± and 
N2 we only have to observe that the bilinear estimates in Corollaries 13.21 and 
13.31 as well as the null form estimate (|3.3() are valid also in the case where both 
u and v solve the homogeneous wave equation, so we can apply the Transfer 
Principle for the H^ e spaces instead of X*'_± h Note also that (|5.17|l must 
be replaced by 

D E U = - £ - 1 7 > (m) + £if " V'^"- 1 ' + £ 4 m " 1) ^ (m " 1) . (6.12) 

The estimate for Nq requires no change. Finally, the estimates for the terms 
involving Aq can be simplified, since we do not care about powers of e here. 
Thus, in N3 we can replace H 1 by L 2 X1 by giving up the e, and then the estimate 
reduces to 1(6.10(1 . Finally, the commutator terms N4 and N5 are replaced by 
a single term, since we do not need to split Aq according to ((5.2(1 . We simply 
expand the commutator and proceed as in the estimate for N$, reducing to 
[((TTUjl and l(6.11|) . This concludes the proof of Theorem IO 
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7 Uniform H 1 bounds and long time existence 

We shall prove: 

Theorem 7.1. Consider the solution (ij) e ,A e ^) of \l-4\l , \l-6)) from Theorem 
\6.1\ existing up to a time T £ > and belonging to the space \6.1\) over this time 
interval. There exist 

(i) a time T* > depending only on sup e>0 H^oIIl 2 ? 

(ii) constants C,M,Eq > independent of E, 

such that if 

XI + Y E <B for all e (7.1) 

then 

X^ + YS<CB for £ < 1 + g and < T < min(T*, T e ). (7.2) 

Moreover, there is a polynomial P with P(0) = 0, independent of e, such that 

Z' T <CZ £ +eP(X^ + Y '), (7.3) 

for T, e as in (J7~2")l 

We claim that this result, together with Theorem 16. II implies Theorem ll.il 
To see this, first observe that the bound in l|7.2|l implies, on account of 1)4. 

W(T)\\ H1 +e K {W{T)\\^ +e\\d t A^T)\\ L2 } < C'B 

for some constant C > C independent of e. Thus Theorem 16.11 implies T e > T* 
for e as in (|7.2|l . In view of the conservation of charge (|1.3f) for the Dirac 
equation, we can iterate this argument any number of times, obtaining 

T E >NT* and X% T , + Yf, T . < (C') N B for s < ^ - {c f° MNgM 
for all N € N. This proves Theorem H~T1 



We shall prove Theorem 17.11 using the iteration scheme from Sect. El In order 
to simplify the notation we drop the superscript e on the fields ip, etc. as 
well as on the XYZ- norms in Definitions 15.11 and 15.21 in the remainder of this 
section, and introduce instead a superscript (m) to denote the m-th iterate of 
a field. We denote by X^ etc. the norms in Definition ! 5 . 1 1 wit h the respective 
fields replaced by their m-th iterate. Then we have: 



30 



Proposition 7.2. There exist C, 7, 6 > and a polynomial P with P(0) = 0, 
all independent of e, such that the estimates 



X 



(rn+l) < + £lp (m) 



T ' 



y^ m+1) < CY + CT A 



z, 



(m) 



(m) 



(m) 



Z^ m+1) < CZ + CT d 



,(m) 



eP 



T ' 

(m) 



T > 



/io/d /or T < 1 and m > — 1, where 



P 



(m) 



fp(x + r ) 

\p(4' 



(m-1) 



(m) 



K. 



(m-l)N 



/or m = — 1, 
for m > 0. 



/n /acf, i/iese estimates hold for ( recall 1^5.1)) ) 

7 < A + 1 - 29. 
The proof is deferred to the end of this section. 



(7.4a) 
(7.4b) 

(7.4c) 
(7.5) 
(7.6) 



Corollary 7.3. There exist C, S > and a polynomial Q, all independent of e, 
such that if 7 > is sufficiently small depending on A, and T, e > are iafcen 
so small that 



2CT s [2C\\tJj q \\ L 2 +1] 2 < 1, 2e^ 2 Q(X + F ) < 1, (7.7) 



then 



<CX + e^ 2 (X + y ), 
Y^ m) <2CTo + £ 7/2 (*o + y ), 



,(m) 



<2Cz + £ 1 -^ 2 (A + r ), 



(7.8a) 
(7.8b) 
(7.8c) 



/or m > . 



Proof. This is a simple induction. Since P(0) = in Proposition 17. 21 there is a 
polynomial Q(r) such that 



P(4[C* + l]r) < rQ{r) for r > 0. 



Then 



P 



(m) 



(7.9) 



holds for m = — 1, in view of the definition (|7.5|l . Hence (|7.8|) for m = follows 
from 17.4a|l - H7.4c|) and the fact that the iterates at m = —1 all vanish. Now 
assume <|7.8ll holds for < m < tuq. Then (|7.9|l holds for such m, and using 
{77) 1 and (f7~4} we obtain lf?~H|) for m = m + 1. □ 
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We are now in a position to prove Theorem 17. II Indeed, from the proof of 
Theorem 16. II we know that the iterates (jf^f* converge in the y-norms. We can 
therefore pass to the limit m — > oo in Corollarv l7.3l Thus, from l|7.8a|) . (|7.8b(l 
we get (|7.2|) . and from (|7.8c|l we get 

Z T <2CZ + 1. (7.10) 

Substituting the latter into the second term in the r.h.s. of (|7.4cll in the limit 
m — ► oo, we then obtain (|7.3|l . This proves Theorem 17. II 

Proof of Proposition \73% By the estimates in Sect. 15.21 

u e± it/^ R (m)U < £ (l/2-A)- p M (7n) 

\\R {m) \\ LHST) <eP^\ (7.12) 

the only difference being that l|5.17|) must be replaced by Q6.12p . Then l|7.4a|l 
follows from the equation (|6.2c|) by applying Lemmas 14 . 21 and l4 . 31 the embedding 
H4.7fl . and the estimates proved in Sect. 15.31 However, instead of the estimate 
l|2.6[) . which was used to prove l|5.23|l . we use the analogous estimate for the 
iterates: 

||n°_^ m) IL, + ||n°^"°|L s < e\\^\\ H1 + £ 2 ||a(™- 1 )||^ 1 ||^™- 1 )|| h1 . 

Next, applying Lemma |4~T1 to the equation (|6.4|l and using the embedding (|4.7|l 
and the estimate (|7.11(l . as well as Ij2.5|l at t = 0, we reduce l|7.4bj) to proving 

II Am) A m )\\ < rfrvWiVW _i_ A~ iv(" l )l 3 

\\ A <t>± H 1 {St) ~ I T J Y T + £ l Y T J • 

But this follows from Leibniz' rule, Holder's inequality, (|5.1U|I and l|5.3|l - (15.61) . 
in view of l|6.2b|l . The factor T s comes from applying Holder's inequality in 
time. Finally, consider l|7.4c|l . Apply Proposition 13 . 51 to (|6.4|l and use (|2.5(l at 

t = to get 

7 (m+l) < y i J-AyWyfm) i V^ll j(™)j,(™)|| i llfi>( m )ll 

T ^ Z + £ T Y T +2^11^0 9± \\ L \ + Lr(S T ) + II-" \\LlLl(S T y 

± 

The last term is covered by (|7.12ll . On account of l|5.5|) and 15.611 , 

II ,(m),(m)|| < ^<5r 7 (™)l 2 || i(™)|| , JyHl 3 

||A) <p± || L i+ Lr(ST) ^ J l z t J II 0± IUj»lj(s t ) + e L y T J ■ 

Then l|7.4c(l follows, in view of 

||j,( m )|| < y( m ) , c \r( m ) 

\\P± \\LfLl{S T ) ~ T -t- fcl T ' 

which holds by (|4.5|) . □ 
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8 Higher order bounds 

Here we prove bounds for higher order derivatives. For m = 0, 1, 2, . . . , set (cf. 
Definitions 15.11 and I5.2fl 



• 




[m] 


— £A S|a|<m 




• 


Yf\ 


m] 


\a\ <rn ^-^± 


\\d%4>%\\ Ki±h ^ {ST 


• 




[m] 


= e A (l|Va§|| H , 


n +e||afll«»), 


• 


Y<?\ 


m] 


= U £ o\\h^- 





The local well-posedness of DM in these norms for m = was established in 
Sect. El and a standard argument shows that higher regularity persists, i.e. if 
X^[m] + Yq[iti] is finite for some m > 1, then Xj,[m] + Yj.[m] is also finite in 
the interval of existence < T < T e . Here we concentrate on proving bounds 
which are uniform in s. Thus, we shall prove: 

Proposition 8.1. // 

A e [m] +Y Q E [m] = 0(1) (8.1) 

then 

X^[m]+Yf[m]=0(l) (8.2) 
for < T < T £ , where T e is the existence time from Theorem \l.l\ 

We claim there exist C, S, 7 > and polynomials Q, P m — all independent of 
e — such that for < T < 1 and m > 1, 

Xf,[m] < CX^[m] + e 1 Q {X e + F £ ) ■ {Xf,[m] + Yf[m}} (8.3) 

+ P m (X e T [m - 1] + Yf[m - 1]), 
F|[m] < C{Y e [m] +eX e Q [m]Y^[m]} + CT s {l + 2CU \\ L2 } 2 Y^[m] (8.4) 

+ e 7 Q (A £ + F £ ) • {Xf,[m] + F#[m]} 

+ P m (Xf,[m - 1] + F|[m - 1]). 

Granting this for the moment, let us prove Proposition 18 . II by induction on m. 
The case m = of l|8.2f> was proved in Sect.[7| Adding up the inequalities (|8.3|l 
and H8.4f> . we see that if (|8.2|l holds for m — 1, then it also holds for m, provided 
T, e > are so small that (|7.7|l holds. Arguing as in the paragraph following 
Theorem 17. II we iterate this argument to cover the full time interval [0,T £ ]. 

To prove the claim, we apply J2\ a \<m ®% t° the system, and imitate the proof 
of the estimates in Proposition l7 . 21 for m = 0. We single out the top order terms 
where m derivatives fall on one of the fields A, ip or <fr±; these are estimated 
exactly like in the case m — 0. All other terms are lumped together and yield 
the term 

P m (Xf, [m - 1] + F# [m - 1]) . 
We skip the straightforward but tedious details of this argument. 
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9 Estimates for the small component 

In this section we prove that if the "positron part" Ilii/) 6 is small initially, then 
it stays small uniformly in every finite time interval, where "small" means either 
0{e) or 0{e 2 ). Here is the precise result: 

Proposition 9.1. (i) Assume i|8.1[l holds for some m > 0. Then if 

ll n -^lff- = °( £ ) 
holds at time t = 0, it also holds uniformly in every finite time interval, 
(ii) Now replace (|8.1|l by the stronger condition 

11^11^=0(1), ||Va§|| Hm + £ ||af|| ffm =0(1), (9.2) 
as e — > 0. Then if 

llni^H^-! =0(e 2 ) (9.3) 
holds at time t — 0, it also holds uniformly in every finite time interval. 

Let us interpret this result in terms of the lower component of e %t l £ ip £ , 
as in (|1.17|l . We claim that (|9.1|l is equivalent to 

\m Hm =0(e) (9.4) 

while (|9.3|) is equivalent to l|9.4|l and 

||%f || Hm _x = 0(1). (9.5) 

The equivalence of l|9.1|l and (|9.4|l follows from (|1.26[1 . since 

\W\\ Hm+1 =0{\) (9.6) 

on account of Proposition 18.11 To prove the rest of the claim, note that by 
lfl~2Tjl . (|9~3|) is equivalent to 

^d 3 rf = O(e), if +is^a j djX e = 0{s 2 ) in H m ~ x (9.7) 

where \ e i s the upper component of e**/ £ ij) £ , as in (|1.17|l . But by the second 
equation in H1.19() . 

ie 2 dtv e = n £ + ie^a 1 d 3 x £ + 0{e 2 ) in H m ~ l (9.8) 

where we used the fact, proved below, that if l|9.1|l and l|9.2|) hold initially, then 

\\VA £ \\ Hm +s\\d t A £ \\ Hm =0(l) (9.9) 
uniformly in every finite time interval. 
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Proof of Provosition WIY i). In view of 12.211 - (|2.4[) . we can replace Ii £ _ip e by <fi- 
in (|9.1|l . and by (|4.5() it suffices to consider the low frequency part (4>-)iow Set 
(cf. Definition ED 

Zj*[m] = 2_j ll 5 "(^-) low IL?L6nL~L2 (ST) ■ 

a|<m 

Then recalling Proposition 18.11 and using induction on m, it suffices to prove 
that there exist constants C, 5 > and polynomials P m — all independent of 
s— such that for < T < 1 and m > 0, 

Zr[m] < C \\<j>.{t = 0)|| ffm + CT S {C ||^|| i2 + l} 2 Z e T [m] 

+ [e + Z e T \m - 1]} P m {X e T [m] + Y£[m]), 

where by convention Zj\—1] = 0. But this estimate follows by a straightforward 
modification of the proof of the estimate for the Z-norm in Proposition 17.21 
taking into account the bound l|7.10|l . 

Let us now prove l|9.9[) . assuming it holds initially. In view of Lemma |4.2I 
this reduces to proving || J e || i 2 ffm ( Sr ) = 0(1). Split J = J' + J" as in Sect. 15.31 
To estimate J", we proceed as in the proof of (|5.22|l . taking into account the 
higher order bound (|9.6(l . On the other hand, since J' has vector components 
| Re ( a 3 (x e )iow, (?f )iow ), we have 

||<9"J'|| L 2 < - || 9 x(X £ )low|| i | iso ll^(77 £ )low|| irL 2 , 

/3+7=q 

and the r.h.s. is 0(1) for |a| < to on account of l|9.4|) . 14.411 and (|8.2|l . □ 

Proof of Provosition WW ii) . Here we break with our earlier notation, writing 

Then from the Dirac equation, 

id t <t>% - ^r^ £ + + {A e <p e ) = o, (9.io) 

idtf. + ^1<PL + ni {A e ^) = 0, (9.11) 
where A £ = A £ a j + A%. Thus, 

01 = (A £ + l)" 1 e 2 {-id t (jf_ - ni (A E (f> £ )} , (9.12) 
so we reduce ()9.3|) to proving 

11^11^-1=0(1), (9.13) 

1^1(^)11^=0(1). (9.14) 
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The latter follows readily from (|9.9|) and l|9. 6[) . since III is uniformly bounded. 
For later use we also note that (|9.10|) implies 

||ft0+||fl»-i=O(l), (9-15) 

since the symbol of x i J 1 is bounded by |£| 2 . 

To prove l|9.13[l we proceed as in ^ Sect. 4], Consider first the case m = 1. 
Take a time derivative of <|fj . 1 1 11 - then take the imaginary part of its inner 
product with dt<j> E _ and integrate in x. Making use of the self-adjointness of 
X E , IT and A E , and the fact that (IT ) 2 = III , we then obtain 



< {\\VA E \\ Hl \\d t cf> E + \\ L2 + \\Vd t A E \\ L2 W\\ H1 ) \\d« 



Thus, dividing by ||^t0— £2 an d integrating in time, and using the fact that 
H9.13[l holds at time t — (this follows from l|9.11|l and the initial assumptions), 
we reduce <|9.13[1 for m — 1 to proving that the terms inside the parentheses in 
the last inequality above are all O(l) locally uniformly in time. But this follows 
from the bounds (|9.6|) . I|9.9|) and (|9.15(l . Here we use also the fact that dtA^ 
enjoys the same bounds as V^4q, m view of the equation 

Ad t A £ = - div J £ 

which follows from ljl.4b|) and the conservation law dtp E + div J e = 0. 

Following we now proceed by induction on m, starting at m = 1. Thus, 
we apply d t d", where |a| < m — 1, to the equation (|9.11() . and we take the 
imaginary part of its inner product with dtd"4> E _ and integrate in x. Then by a 
straightforward modification of the argument for m = 1, we reduce (I9.13|) to the 
O(l) bounds H9.6(l . I|9.9|) and Ij9.15|l . as well as (|9.13(l at the previous induction 
step. We omit the details. □ 



10 Nonrelativistic limit 

We first prove Theorem 1 1.21 then we discuss the modifications needed to prove 
Theorem 11.31 

Proof of fl.lcia)) . This can be restated: 

e l */ £2 n° r , e- lt / £2 U°_r - v- = ( ° ) in H 1 as e -> (10.1) 

locally uniformly in time. We claim it suffices to prove 

e lt/£2 Il E + ip £ -> ( V +) , e- lt/E2 IL E _ip £ -> v_ = ( J in H 1 as e -> 0. (10.2) 
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To prove the claim, write H £ ± ip £ — H ±ip £ ± r £ . By the orthogonality between 
II+ and n°_ , we get 



„e 



But if (|1U.2I) holds, then the right hand sides converge to zero in if 1 . Thus 
r £ = o(l) in H 1 and we have proved that (|1U.2|) implies l|lU.lfl . In the remainder 
of the proof we skip the superscript e on the fields, to simplify the notation. 
Using (|2.2() and (|2.4|) we reduce H10.2|l to proving 







in H 1 as e 



0. 



(10.3) 



uniformly in any given time interval [0, T]. By (the proof of ) Theorcm ll.il the 
solution exists in this time interval for all sufficiently small e > 0, and 



X§. + Yf. = 0(1), 



\e ±lt/£ R e \\ x i,e-i ,„ > =o(l) (10.4) 



as e — > 0. Note that i|10.3fl holds at time t = 0, by Lemma 11.61 Thus, it 
suffices to prove that there exist K, 5 > 0, depending on T and Xf, + Yj., but 
independent of e, such that for every time interval I = [to, t\] C [0, T], 



f(I)<Kf({t })+K\I\ d f(I)+o(l) 



as e — > 0, where 

m-- 







(10.5) 



(10.6) 



W.l.o.g. we assume I = [0, T], and we only estimate the first term in (|10.6|) . 
Write 



= (7 E (i)0+ + / (7 £ (i - s) ds, 



v+(t) = S(*H + 



S^f — s) [(uv+)(s)] ds, 



where </>q , Vq are the data of <fi+, v+ and U £ {t), S(t) are given by (|4.1(l . Thus 



0+W 



(t) = c/ E (t) 



+ [C/ £ (t)-5(t)]( 1 ' 



fu £ (t-s) 
Jo 







(,s) + L £ + <M S ) 



(10.7) 



+ / [S(t-s)~U £ {t-s)] 
Jo 

= h + h + h + h- 







(s) (is 
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Clearly, 



\M\L^m(S T ) ~ 



if 1 



As in El Sect. 5], 

W^Wlfh^St) = for J = 2 ' 4 > 

using the dominated convergence theorem and the fact that 



IVuI 



LI 



< 



i + ll^-llffi < 00 



(10.8) 



uniformly in every finite time interval. It remains to consider I3. By Lemma 
EUand the embeddings and ^Jj , 



< 







.4 



0<P+ 



,±it/s 2 R s\ 



A 



(10.9) 



The second term on the r.h.s. is o(l) by l|10.4(l . and the first term is bounded 

by 



T 



1/2 



L^H 1 (St) 



+ \\(U- ^o)^+|lz,oc H i(s T ) 



But using Leibniz' rule, Holder's inequality and Sobolev embedding, it is easy 
to see that the terms inside the parentheses are dominated by K /(J), where K 
depends on the size of Xj. + Yj, and (|lt).8(l . □ 



Proof of f 1.1 3b]) and \1. 13c]) . Using Sobolev embedding we reduce l|1.13b(l to 
H1.13c|) . To prove the latter, observe that (|10.3|) implies 



X-\ 



v+, X- 



0, rj + — > 0, r\- — > v_ in H 1 as e — > 0, 



(10.10) 



locally uniformly in time. Thus (|1.13cfl follows immediately from (|2.15ll using 
Holder's inequality and Sobolev embedding. □ 

Proof of $1.15]) . Multiply l|2.16[l by a C 1 compactly supported test function 
G(t, x) and integrate in t, x. W.l.o.g. assume G is real- valued. The integrals 
corresponding to the last two terms in r.h.s. 1)2. 16[) are 0(e) in absolute value. 
To see this, integrate by parts in time and use 



fgdx 



< \\9\\m 



and the bound, locally uniform in time, 



0(1). 



(10.11) 



(10.12) 



The latter is easily reduced to the uniform bounds for Xj, + Yf,, using Lemma 
12.11 Sobolev embedding and Holder's inequality. 
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Next, fix 1 < j < 3 and consider 

I± := - Re J (a^ x±,V±) G dtdx. 

In view of fTITTrHl and 

7_ = ^Rc J (o 3 X ~,v^)Gdtdx + o{\). (10.13) 
By (|22J),(E3I) and dTTT^Il . 

~ = e - lt / £2 in°n e _^ + o(£ 2 - A ) in l 2 x (io.i4) 

locally uniformly in time. But by (|1.28|) . 

1^. i M - (•**<) + ± (!-[»• 
In view of l|1.13a|l and the bound 11.29J1 . it follows that 

= ^[AT 1 (^) + i (1 - [AT 1 ) ( 6 " 2 t V ) + 0(1) (10.15) 

in L 2 . Moreover, by dominated convergence, 

[\ £ ]~ 1 a k d k v- = cr k d k v- + o(l) in 7T 1 . (10.16) 

Using (|10.13|) (|1U.1()[) and either Holder's inequality or (|10.11f> . we conclude 
that 



where 



I_=Rc J i(cr j a k d k v^,v^)Gdtdx + I'_ + o(l) 
71 = ^Re /" e- 2 " /e2 (^'(l-[A e ]- 1 )« +) i;_>Gdida;. 



But the latter is 0(e) in absolute value (integrate by parts in time and use the 
analogue of (|10.12fl for v±). Using i|1.2|) we finally conclude that 



I- = J ^-Im(d ] v^,v-)-^ kl d k (a l V-,v_} 



Gdtdx + o(l). 



A similar calculation can be done for 7+, and this proves (|1.15(l . □ 

Next, we prove Theorem 11.31 By hypothesis, (|9.4(l . or equivalently 19.1(1 . 
holds initially and therefore also uniformly in every finite time interval, by 
Proposition 19. II Next observe that since (|1.16l) holds initially, we have 

B4V £ = V4 + O(e) in H 1 
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locally uniformly in time. In fact, this follows from H2.2f> - H2.4f> . since (|8.2[l holds 
with m = 1. Wc conclude that it suffices to prove l|1.16fl with i/> £ replaced by 
if)^_. We proceed as in the proof of Theorem ll.2l but now the remainder term in 
<|lU.5fl must be improved from o(l) to 0(e), and /(/) is given by the first term 
in r.h.s. l|10.(j[l . Again we reduce to estimating the terms as given by 

(JTUTJl. 

The term Ii is estimated exactly as before, but is now 0(e) since (|1.16|) is 
assumed to hold initially. Using the fact that U e (t) — S{t) — e 2 lZ\(t), where 
7£f(t) is bounded from H s+A — > H s uniformly in e and < t < T, and the 
assumption that the initial datum of w + is in H , we find that 

For I3 we use again l[l().9jl . but now the last term is 0(e), as follows from the 
proof of Ij5.12|) taking into account the fact that (|8.2|1 holds with m = 1. The 
first term in r.h.s. l|lU.9|) is estimated exactly as before. This proves and 
then it follows immediately that l|1.13b(l and i|1.13c|) are also improved to O(e). 



11 Semi-nonrelativistic limit 

Here we prove Theorem 1 1.71 The initial assumptions (i), (ii) imply, as proved in 
Sects. El and H that (JH3 holds with m = 4, while (|nTJ-<EII> hold with m = 2 - 
We write 



:=e 4 * /e V 



is ._ Jt/e 2 ,e 



with ip± defined as in (|2.1|l . Also, we denote by x± the upper component of </>§_. 
Observe that 

=^±+0(e 2 ) in H\ 

in view of (|2.2|) - (|2.4|) and the bounds (|9.6|) and (|9.9I) for m = 2. On account of 
1)9. 3p . we may therefore replace x £ i n Hl-31|) by x+- By Lemma \'Z. II 



id t 



A £ - 1 



where 



Rf = \eiA^ ■ V X E - \eB^ X e + \e 2 (A £ ) 2 



~ (1 - [AT 1 ) e {2*A e • V X £ - BjtrV} + ^[A e ] _1 £ {^V } 



i(l-[A1- 1 )e 2 {(A^) 2 x £ }-^ 2 [A1- 1 



.4 



A e - 1 



(X+ - X- 



Recalling the bound l|1.29|l on the symbol of 1 — [X 6 ] 1 and using the fact that 
(j8.2j) . (|9.4j) and (|9.9|) hold with m = 2, we conclude that 



£ E = -«A E • V X £ - -eB^ X e + - £ 2 (A £ ) 2 X e + 0{e 2 ) in if 1 , 
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locally unformly in time. 

Then, since (|8.2I) holds with m = 4 and ^± = y + £ 2 7?.4, where TZ% is 
bounded from H s+A — > _ff s uniformly in e, we further conclude that 

= 5 (iV + sA^f x% ~ Alx% - \eBy X % + e 2 r s (11.1) 

where r e = 0(1) in H 1 locally unformly in time. Comparing Ull.lfl to the Pauli 
equation (|1.30(l via the energy inequality for the self-adjoint "Pauli operator" , 



one finds that 
as e — > 0, where 



f(I)<f({t a })+K\I\f(I) + 0(s 2 ) 



f{I) \\X+ Xp\\ L *> H i( IxV a} 

for time intervals / = [t , t{\ C [0, T], and where K depends on T but not on e. 
In fact, K depends on the O(l) bounds in (|9.6|) and (|9.9|l . which hold for m = 2 
as we recall. We conclude that /([Q,T]) = 0(e 2 ), and this proves l|1.31fl . 

Observe that i|1.23ll holds in H 1 locally uniformly in time, in view of l|1.20[l 
and the fact that (|fj.4|> . I9.5jl and (j9.9|l hold for m = 2. Substituting (|1.23|) into 



J-= £ - 1 {2Re(aV,^) c2 } fc=1 



2.3 



and using lfT3T|) yields lfT32Jl . 

12 Proofs of the spacetime estimates 

Here we prove Theorem 13. II and Proposition ^. 41 

Proof of Proposition Let Q be a cube with side length ~ \x centered at £ j 
where |£o| ~ A, and let Xq(0 be a smooth cut-off function equal to 1 on Q. For 
example, we can take 

xq(0 : =^— 



(12.1) 



where r\ is a smooth bump function equal to 1 on a neighborhood of the origin. 
Then by the TT* method, we reduce (|3.7|l to the decay estimate 

|tf.Q(t,*)l<( M|t| M i f ° rA - 1/ / £ ' (12.2) 
i e,yv ;i ~ I^Ajip 1 for A> 1/e, v ; 

for the convolution kernel 
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with h e given by (|2.14|) . In view of the scaling identity 

K StQ (t,x) = e^K^qie-h^^x), 

it suffices to prove <|12.2[1 for e = 1. To simplify the notation we write Kq 
instead of K\^q. Thus, 

K Q {t,x)= r I e lrx ^e lta ^ XQ {ruj)r 2 da{u:)dr (12.3) 
Jo Js 2 

where a is surface measure on S 2 and a is given by l|2.18|l . Note that 

a 'W = 7T=f and Q ^)= (1 + ' 2)3/2 - (12-4) 

We split the problem into the following cases: 

(i) A < 1 and \x\ >X\t\, 

(ii) A < 1 and \x\ < X\t\, 
(hi) A > 1 and |x| > |t|, 
(iv) A > 1 and |x| < 

Rewrite (|T2~3t as if Q (i, cc) = J °° e l * a Ma(r, x)r 2 dr where 

a(r,x):= f X qM^H- 
We shall need the following: 

Lemma 12.1. |a(r, a;)| < (r\x\)~ 1 \i(r), where xi is the characteristic function 
of an interval I of length ~ /i and centered at a distance ~ A from the origin. 

Proof. The statement about the r-support of a(r,x) is obvious, and the decay 
statement follows from the fact that 



e lx ' ul 1 (uj) da(oj) 



<i/N 



for all smooth functions 7 such that I7I < 1. But this fact is easily proved by 
passing to spherical coordinates and rescaling. □ 

Thus 

\K Q (t,x)\<J(r/\x\)dr~n\/\x\, 

and this covers the cases (0) and (JmJ) above. 

To handle the remaining cases we write Ijl2.3|l as Kq(t, x) — J s2 b(uj) da(oj), 
where 

d 



b(u) 



dr 



^i(ta.(r)-\-rx ■ uj) 



? X ? Mr , dr. 
iyta'yr) + x ■ w) 
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Integrate by parts and write 



d_ 

dr 



2 



i(ta'(r) + x ■ uj) 



X Q (rui)r 2 ta"(r) £[x Q (ruj)j 



i(ta'(r) + x ■ uj) 2 i(ta'(r) + x ■ w) 



Correspondingly we split b = b\ + 62. Observe that the r-support of XQ( ruJ ) is 
contained in an interval I of length ~ /i and centered at a distance A from the 
origin, while the w-support is contained in a set given by 

Z( w ,w )< M /A (12.5) 

for some wo G S 2 . Moreover, in view of (|12.1|1 we have 

\£x Q (ru>)\<l/fi. (12.6) 

Now consider case Jrv}. Then on account of (|12.4|) we have a'(r) ~ 1 and 
a"(r) ~ A -3 for r S I, so \ta'(r) + x-uj\>\t\. Thus 



|6i(o;)|<(l/A|t|) J dr<[i/\\t\, 
which is more than good enough. Next, using l|12.(jfl we have 

I&2MI < (V \t\) [{r + r 2 /v) dr < (A/i + A 2 )/ |t| < A 2 / |t| 



But integrating this over the region (|12.5|) on S 2 gives us a bound fi 2 / \t\, which 
again is more than good enough. 

Finally, consider case JuJ. Then a!{r) ~ A and a"(r) ~ 1 for r £ /, so 
|ta'(r) + x-u>\>\\t\. Thus 



|6iHI<(V|t|) J^rZu/m 

and 

IfcMI < (l/A |t|) jf (r + r 2 /y) dr < + A)/ |t| < A/ |t| . 

Taking into account (|12.5|) we thus get the desired bound, and this concludes 
the proof of Proposition [jO] □ 

Proof of Theorem \3 . ltivfy . If /1 ~ A, this reduces to part (fTTTfl of the theorem, so 
we may assume /1 <C A (and A > 1/e) . But then by an orthogonality argument 
(see, e.g., the proof of the analogous estimate in Theorem 12.1 of |BJ) we reduce 
to proving 

in the case where the Fourier transforms of /, g are supported in (diametrically 
opposite) cubes with side length ~ /1 and at distance ~ A from the origin. 
But this follows from Holder's inequality and the estimates l|3.5|) and (|3.7() with 
(g,r) = (4,4). □ 
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Proof of Theorem Vd '. Jlf z)) . If \l ~ A, this reduces to part JmJ) of the theorem, so 
we may assume /i<A< 1/e. By orthogonality, we reduce to proving 

IIHI^^Vll/MMI^ (12.7) 

in the case where /, 17 are supported in opposite cubes Q, —Q with side length 
~ \x and at distance ~ A from the origin. By rescaling t — ► i/e we further reduce 
to proving l|12.7fl without the e 1 / 2 in the right hand side, and with u, v given 

by 

HtTiO = e« lel /(£), MtTiO = e^^giO- (12.8) 
Here a is given by (|2.18J) . Then by a standard Cauchy-Schwarz argument, see 
e.g. Sect. 3.4], we finally reduce to proving that 

X{r,-v£Q}n{rr.ri-teQ}(v)ti(T ~ \v\ ± - »7|)) dij < m 2 (12.9) 
where 

k(p) := e - 1 a(ep). (12.10) 

(Here and in what follows we use the notation %a for the characteristic function 
of a set A.) Then in view of i|12.4[l there is an absolute constant Co such that 

\k'(p)\ < Co < 1 for all p < 1/e, < e < 1. (12.11) 

Denote by I±(r, £) the integral in l|12.9|l . In polar coordinates 77 = rw, r > 0, 
w € .S 2 , we have -Z±(t,£) = J s2 a±("r, £; w) da(u>), where 

/•oo 

a±(T,£;w):= / X{r:™eQ}n{r : ™-CGQ}(r)<5(r-r±/c(|^-rc i ;|))r 2 dr. 
Jo 

Observe that the w-support of a is contained in a set given by Ijl2.5|l . so it 
suffices to prove that a± < A 2 . Observe also that in the integral defining a±, 
the variable r is restricted to an interval / of length ~ fi and centered at a 
distance A from the origin. 

We shall use the following fact: If / : R — > R is differentiable with |/'(r)| > 0, 
and / has a zero at ro, then 

Take 

/(r):=T-r±fc(|£-w|), (12.13) 
for fixed r, £, cj. Then for r such that - £ £ Q, 

\f'(r)\ = 1 T - H ^r^r > 1 - Co > 1, (12.14) 

where we used (|12.11|) and the assumption A < 1/e. On account of 112.1211 and 
H12.14|l . we then get a± < A 2 as desired. This concludes the proof of part (0 of 
Theorem O □ 
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Proof of Theorem Iff, ilf m|) . This reduces to proving 

X{r,:\<n\^^}n{r,:\i-n\^X}{v)^{T ~ M ± K\£ ~ v\)) d V ^ [min(^. A)] 2 , (12.15) 



for k defined by l|12.10|) . Let us denote the above integral by I±(t, £). Passing 
to polar coordinates we have I±(t, £) = J g2 a±(r, dcr(w), where now 

/>0O 

a±(r,£;w):= / X{r:r~ M }n{r:|™-£|~A}M<5(T - r ± - rw|))r 2 dr. 
Jo 

We split into the cases 

(a) A < 1/e, 

(b) A > 1/e. 

Case Then in view of (|12.12[) and (|12.14|l with f(r) given by (|12.13|) . we 

have a± < /x 2 . Now integrate over 5 12 , taking into account the fact that on the 
support of a± , 

=Av,t) Z A/a* if A^>A. (12.16) 

Case By rotational symmetry we may assume £ = (|£|,0,0). Now 

parametrize the sphere S 2 by 

(y,0)»u= (y^l-y 2 ra(0)) , fl(0) - (cos 0, sin 0). 

Then surface measure da(ui) on S 12 becomes dydd. Again we use (|12.12fl with 
f(r) given by 112.13(1 . Observe that / depends implicitly on y but not on 8. 
Denote by A = A(t, £) the set of y € (—1, 1) such that f(r) given by l(12.13[) has a 
zero ro = ro(y) > 0. Since |/'(r)| > 0, the implicit function theorem guarantees 
that A is open and ro : A — > (0, oo) is a smooth function. Differentiating 
f( r o{y)) = gives 

= f'(r Q y (y) T m ~ r u\)-f^- (12.17) 

where we used £ • c^w = £i = |£| and u> ■ d y u! = 0. 

Let us suppress the subscript and write r(y) instead of ro(y) from now on. 
Solving (|12.17|) for r'(y) and using the fact that f'(r) < 0, we see that dr/dy 
is either strictly negative or strictly positive, depending on whether we have 
the + sign or the — sign in l|12.15fl . The function r(y) is therefore a change of 
variables. 

With this information in hand, we solve (|12.17fl for f'(r) and substitute into 
<|12.12|) . thus arriving at the identity 



F(r,)6(T-\ V \±k(\£-r ) \))dr l = I I F(ru), '' K '' 4 



m>{\Z-ru:\) 



Or 
dy 



dydO. 
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Changing variables y — ► r finally gives 

F(t])S(t - \r)\ ± fc(|£ - dr) = [ [ F(ru) T |f~ — drdO, (12.18) 



where w is now a function of r and 0. We apply this with 

Pi 7 !) : = X{? 7 :|'7l~A*}n{ij:|e-T7|-A} (^)- 

Since A 3> 1/e, we see from (112.4(1 that — ru>|) ~ 1, whence 

- (12.19) 

We now split into the subcases 
(bl) /i«A, 
(b2) M ~ A, 
(b3) (i > A. 

Case In this case we can prove the estimate in Theorem 13. 1 Ifiiip) directly, 

by applying Holder's inequality followed by the linear Strichartz estimate (|3.6fl 
with (q, r) = (4,4). (This works because we are at high frequency, i.e. 3> 1/e.) 

Case (bl). Then |£| ~ A, so the desired estimate <|12.15ll follows readily from 

ima and trrm . 

Case f&3). Then |£| - ^, so (tT2~lH and iT2~T£l) imply 



J±(t,£)<A / / X{T:T~n}r\{r:\i-Te>\~\}{r)drdO. 

Recall that uj is now a function of (r, 9). However, |£ — ru>\ is independent of 9, 
so by a slight abuse of notation we will simply write lu — uj(r) and integrate out 
9, leaving us with 



A / X{r:r~n}n{r:\£-rLj(r)\~\}{r) dr. 

Clearly it suffices to prove that the support of the integrand is contained in an 
interval of length ~ A. Let us assume there is no such interval, and obtain a 
contradiction. Fix a point ro in the support, and write 

r — ro + ft 

for a general point r in the support. In view of our assumption, k varies on a 
scale 3> A. Thus, if we can show that 

|£ - ruj(r)\ 2 = a + K 2 + 0{\n + A 2 ), (12.20) 



4G 



for some constant a, it follows that |£ — rui(r)\ also varies on a scale ^ A, and 
we have the contradiction we seek, since |£ — ruj{r)\ ~ A on the support. 
To prove H12.2l)[l . write 

|e-rc| 2 = |e| 2 + (r 2 -2r|e|) + 2r(l~ Wl )iei- 

On account of (|12.16|) we have 1 — uj\ < (A//1) 2 , so the last term on the right 
hand side is 0(A 2 ). For the second term we calculate 

r 2 - 2r \i\ = (rg - 2r |£|) + 2(r - |£|)k + * 2 . 

But 

|r -|el| < kow(ro) - £| ~ A, 
so we conclude that (|12.2U|I holds. This ends the proof of Theorem 13. II □ 

13 Proof of Theorem [O 

As remarked, by a standard procedure this reduces to some well-known bilinear 
estimates for the homogeneous wave equation. The first observation is that 
by rescaling x — > ex we can reduce to the case e = 1. Thus we suppress the 
subscript on H s ' e etc. from now on. 

Some notation: For s € K, let D s , D s + and D s _ be the Fourier multipliers 

(D s u)~ = \£\ s u, (Dlu)~={\T\ + \t\) S u, (Dluy=\\r\-\(;\\ s u. 

The notation u ^ v means \u\ < v. We are concerned with bilinear operators 
B(u, v) of the form 

[B(u, v)r(r, = J b(r - A, £ - 77; A, t])u{t - A, £ - ry)w(A, 77) dA dry, 

where 6(t, £; A, 77) is the symbol of B. The symbols of the null forms Qo, Qij 
and Qoj are, respectively, 

9o(t,£; A, 77) = tA - £ • 77, (13.1a) 
<7-y (r, f ; A, 77) = -^77j + Cj77i, (13.1b) 
?0i(r,C; A, 77) = -TT/j + A£,-. (13.1c) 

Since we rely on estimates for the absolute values of these symbols, and since 
all norms involved only depend on the absolute value of the Fourier transform, 
we may assume u, v > henceforth. 

For s £ R, let R s be the bilinear operator with symbol r s , where 



Kt,& A, 77) 



\Z\ + \r)\-\Z + v\ if ^A>0, 
\Z + V\-\\Z\-\V\\ if ^A<0. 
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We shall need the estimate, for 9 > 1/2, 

\\R 1/2 {u,v)\\ L2 <\\u\\ H0 , s \\v\\ Hm , 9: (13.2) 

which derives from an estimate for the homogeneous wave equation via the 
Transfer Principle; see JS] for the details. We also need 

R s (u,v) ^ D s _(uv) + {D s _u)v + uD s _v (s > 0). (13.3) 

This follows easily from the triangle inequality, if one keeps track of the signs of 
r and A as in the proof of the following lemma, which is more or less standard. 

Lemma 13.1. The following estimates hold: 

Q lJ (u,v) < R 1/2 {Du,D 1/2 v) + R 1/2 (D 1/2 u,v) (13.4a) 

* RW (Du, D^ + v) + (D^ + u, Dv), (13.4b) 

Qoj(u, v) £ [r./t.s. (|13.4a|l ] + Du ■ D_v + D_u ■ Dv, (13.4c) 
Qo(u, v) ^ [r./t.s. (|13.4all ] + D+u ■ D_v + D_it • D+v. (13.4d) 

Proof. All these statements reduce to estimates on the absolute values of the 
symbols (113. Q . First, by Lemma 13.2] we have 

|g«(r,f;A,»j)| < |exr?| < \t\ 1/2 \v\ 1/2 \t + v\ 1/2 [r(r, £ A, V )] 1 ' 2 , 

where r is the symbol of R as defined above. Then (|13.4a|) and (|13.4b|) follow, 
in view of the fact that 

r(T,£;A,77) <2min(|£|,M)- (13.5) 

To prove l|13.4cl) . write 

q j{T,Z;\,T]) = (ei |£| - r)r/j + (A - e 2 -ei(|£|?7j - £1^2 

where ei and e 2 are the signs of r and A, respectively. That is, t\T = |r| and 
e 2 A = |A|. Now take absolute values and use the fact (see jSJ Lemma 13.2]) that 

| \t\ Vj ± \r,\ & | < |C| 1/2 \v\ 1/2 + M) 1/2 [r(r, £; A, ri)} 1/2 

holds for all t, £, A, r\. (The sign in the left hand side is independent of the signs 
of r, A.) This proves (|13.4c|l . The proof of (|13.4d|) is similar. Write 

qo(r, £ A, rj) = {r- ej |$|)A + (A - e 2 l^ex |£| + £l e 2 M |£| - £ • r?. 

Then use (see jS] Lemma 13.2]) 

I l€l - ^ - »?| < (lei + l»?|)r(T-, C; A, ii) 
and (JT33|1 . □ 
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Finally, we need the estimate (here si, Sa, 6\, 02 > 0) 

\\ uv \\l 2 ~ \\ u \\h°i-»i IMI.fp2.e2 f° r s 1 +s 2 >| ) 6*1+6*2 > ^- (13.6) 
See Proposition A.l] for the simple proof of this fact. 

We are now ready to prove Theorem 15.31 By interpolation, we reduce to 

\\Q( u i v )\\ L 2 - IMImm IMI-h2,£) ( 13 - 7 ) 

\\Q(u,v)\\ H0 ^ 1/2) - < ||u||^ M ||«|| w( 3/2)+, 1 (13.8) 
where ||w||^i,e in the right hand side can be replaced by ||w||^i,« if Q = Qij- 

Proof of l|13.7|l . First observe that for the last two terms in the right hand sides 
of (|13.4c|l and (|13.4d[) . the estimate reduces to special cases of (|13.6|l . since we 
can always replace D_ by D e _D^T e . Thus, it only remains to prove the estimate 



for the right hand side of l|13.4a|l . but this reduces to 1)13. 2|) . □ 



Proof of i|13.8fl . First consider Qy. Applying l|13.3|l to (|13.4b|l . we reduce to 

IIHL* ^ IMIffo, (1 / 3 )+ , (13.9a) 

Il™lli2 < ||u|| J j(i/ 2) -,(i/2)+ ||w|| ff( i/2)+,i , (13.9b) 
||™|| ff o. ( -i/2)- < \\u\\ L 2 \\v\\ Hl +,i , (13.9c) 
HHI ff o.(-i/2)- < ||u||h(i/2)-,o ||u|| H (i/2)+,i , (13.9d) 
||ww|| ff o, ( -i/2)- < ||M|| ff o, ( i/2)+ IMI H i+, ( i/2)+ , (13. 9e) 

ll UV ll_f/0,(-l/2)- < ||u||jj(l/2)-,(l/2)+ 11^11^(1/3)+, (1/2)+ ■ (13. 9f) 



Via duality and the Transfer Principle, these reduce to the estimates in Corol- 
laries 13.21 and 13.31 which are valid in the case where it, v are both solutions of 
the homogeneous wave equation, as remarked in Sect. [31 

It remains to consider the second and third terms in the right hand sides of 
l|13.4c|) and (|13.4d|l . For the second term we can apply l|13.6fl directly, while for 

the third term we replace D- by D^ 2 ' D^ 2 " 1 , thus reducing to (|13.9d(l . □ 
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